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Preface 


The  purpose  of  this  work  was  to  gather  data  and 
calculate  Franck-Condon  Factors  for  several  of  the  diatomic 
halogen  and  interhalogen  molecules.  Currently  there  is  a 
renewed  interest  in  these  values  with  the  hope  that  they  may 
both  aid  the  spectroscopist  in  understanding  experimental 
observations  and  aid  research  into  new  lasing  systems  by 
providing  early  knowledge  of  the  probabilities  of  transitions 
occurring  between  excited  states  of  these  molecules. 

I  enjoyed  working  on  this  project  and  I  feel  that  I  have 
learned  a  great  deal  through  doing  it.  I  would  especially 
like  to  thank  my  advisor.  Dr.  Ernest  A.  Dorko,  for  taking  the 
time  to  explain  all  of  the  required  background  material  to 
me.  Doing  a  topic  outside  of  my  class  area  meant  that  I  had 
many  new  and  strange  ideas  to  learn.  With  Dr.  Dorko's  help, 

I  gained  a  great  deal  of  understanding  in  a  subject  new  to 
me . 

One  of  the  most  challenging  and  frustrating  aspects  of 
this  thesis  was  learning  the  computer  operating  system  and 
running  the  programs.  I  would  like  to  thank  Major  Tom 
Hopkins  for  his  good-natured  help  and  for  answering  all  of  my 
many  computer  questions.  Finally,  I  would  like  to  thank  my 
husband  and  fellow  classmate,  for  going  through  it  all  with 
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Abstract 


This  study  calculated  Franck-Condon  Factors  for  several 
diatomic  molecules.  For  one  of  the  molecules,  lead  oxide, 
the  molecular  constants  required  to  do  the  calculations  were 
obtained  from  experimental  observations  of  the  compound 
performed  at  AFIT.  Potential  energy  curves  were  generated 
for  eight  excited  states  (a,  b.  A,  B,  C,  C',  D,  E)  and  the 
ground  state  (X).  Franck-Condon  Factors  were  calculated  for 
the  transitions  of  each  excited  state  to  the  ground  state  as 
well  as  the  A  to  a  and  D  to  a  transitions.  Data  obtained 
from  the  open  literature  was  used  for  the  remaining  five 
halogens  or  interhalogens  studied;  molecular  bromine,  bromine 
chloride,  bromine  monofluoride,  iodine  chloride,  iodine 
monofluoride.  The  electronic  states  modelled  and  the 
transitions  between  states  for  which  Franck-Condon  Factors 
were  calculated  varied  with  each  molecule. 

The  potential  energy  curves  were  modelled  using  a 
program  written  by  Dr.  C.  R.  Vidal  of  the  Max-Planck 
Institute  in  Gartach,  Germany.  The  program  uses  a 
Rydberg-Klein-Rees  approximation  to  initially  model  the 
curves.  These  values  are  then  shifted  and  refined  using  an 
Inverted  Perturbation  Approach  routine.  This  is  an  iterative 
procedure  which  increases  the  accuracy  of  the  final 
potential,  especially  in  the  higher  energy  regions.  The 
RKR-IPA  program  also  solves  the  Schrodinger  Wave  Equation  for 
each  vibrational  level  and  creates  a  file  containing  a  list 
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of  wavefunction  values  for  each  level  at  a  specified  number 
of  points  along  the  interval  of  internuclear  separation  of 
interest.  The  Franck-Condon  Factors  are  calculated  using  a 
program  written  by  Capt  Joseph  Pow  and  modified  in  this  work 
to  upgrade  the  approximation  of  the  overlap  integral  from  a 
simple  histogram  approach  to  a  Simpson's  Rule  approximation. 
The  Franck-Condon  Factors  calculated  in  this  work  were  then 
compared  to  values  found  in  the  literature.  Good  agreement 
between  the  most  probable  transitions  (the  Condon  Parabola) 
of  each  set  of  factors  was  found  even  though  the  numerical 
values  of  the  two  sets  differed.  This  shows  that  the  less 
accurate  RKR  approach  is  sufficient  for  many  uses  of 
Franck-Condon  Factors. 


CALCULATION  OF  FRANCK-CONDON  FACTORS 


FOR  DIATOMIC  MOLECULES 

I.  INTRODUCTION: 

Many  properties  of  a  diatomic  molecule  can  be  described 
by  a  set  of  potential  well  functions,  one  for  each  electronic 
state  of  the  molecule.  Transitions  between  two  states  can  be 
described  by  looking  at  a  spectrum  or  the  properties  of  the 
transition  can  be  modelled  by  considering  the  potentials  of 
the  two  states.  When  the  wavef unctions  for  vibrational 
levels  within  the  electronic  states  are  combined  through  the 
Franck-Condon  Factor  the  probability  of  that  transition 
occurring  can  be  predicted.  This  is  especially  useful  to 
those  developing  new  lasing  systems  which  depend  heavily  upon 
optimizing  transitions  from  select  molecular  energy  states. 
Prior  knowledge  of  the  probabilities  of  such  transitions  can 
allow  more  knowledgable  selection  of  the  materials  to  be 
tested . 

The  calculation  of  these  Franck-Condon  Factors  makes  use 
of  molecular  constants  which  are  determined  from  experimental 
observations  made  by  spectroscopists .  These  constants  are 
used  in  the  calculation  of  the  potential  well  for  the 
electronic  state.  This  potential  is  then  used  to  solve  the 
Schroedinger  equation  to  give  the  energy  and  the  wavefunction 
for  each  vibrational  energy  level.  The  two  wavef unctions 


(one  from  each  state)  are  multiplied  together  and  integrated 
over  the  region  where  they  overlap.  The  square  of  this 
integral  is  the  Franck-Condon  Factor.  Spectroscopists  also 
use  the  factors  to  confirm  or  deny  their  transition 
assignments  to  lines  in  an  observed  spectrum. 

This  thesis  effort  will  concentrate  on  the  calculation 
of  Franck-Condon  Factors  for  diatomic  molecules.  In 
particular,  PbO  will  be  investigated  using  molecular 
constants  determined  in  a  recent  thesis  effort  (11).  PbO  was 
previously  investigated  in  a  thesis  by  Capt  Joseph  J.  Pow 
(20).  Results  of  the  two  calculations  are  compared  for  one 
of  the  states  of  lead  oxide.  In  addition,  Franck-Condon 
Factors  will  be  generated  for  interhalogens  of  interest;  IC1, 
IF,  Br  (molecular),  BrCl,  and  BrF.  The  potentials  and 
wavef unctions  were  generated  by  a  program  originally  written 
by  Dr.  C.R.  Vidal  of  the  Max-Planck  Institut  fur  Physik  und 
Astrophysik  in  Gartach,  Germany.  This  program  makes  use  of 
the  Rydberg-Klein-Rees  (RKR)  method  to  approximate  the 
potential.  The  program  also  includes  an  Inverse 
Perturbational  Analysis  (IPA)  to  "fit"  the  potential  curve 
better  at  the  higher  energy  values  (10:  22  Aug).  The  IPA 
routine  is  optional  and  may  be  bypassed  if  desired.  Vidal's 
program  was  written  to  run  on  a  Cray  computer  and  therefore 
had  to  be  adapted  to  run  on  the  CDC  Cyber  computer  available 
to  AFIT  through  the  Aeronautical  Systems  Division.  This  was 
done  by  Capt  Lyle  L.  Rutger  in  his  1983  thesis  (22).  Using 
the  output  from  Rutger's  efforts,  Pow  wrote  a  program  to 


< 


correctly  allign  the  wavef unctions  and  generate  Franck-Condon 
Factors.  The  integration  routine  was  upgraded  to  a  Simpson's 
Rule  approximation  from  a  histogram  approach  for  the  present 


thesis . 


II.  LITERATURE  AND  BACKGROUND  THEORY: 

The  generation  of  the  vibrational  wavef unctions  is 
central  to  calculation  of  Franck-Condon  Factors.  In  order  to 
do  this  a  numerical  solution  to  the  Schrodinger  Wave  Equation 
must  first  be  found.  This  in  turn  requires  knowledge  about 
the  form  of  the  potential,  U(x).  The  potential  model  used 
for  these  calculations  is  a  Rydberg-Klein-Rees  (RKR) 
approximation  with  an  Inverted  Perturbation  Analysis  (IPA) 
used  to  increase  the  accuracy  of  the  model  at  higher  energy 
states.  The  RKR  approach  is  used  when  the  molecular 
constants  for  each  energy  state  are  known  and  the  form  of  the 
potential  is  desired.  Values  for  the  constants  can  be 
obtained  from  any  source  and  are  expressed  in  units  of 
inverse  centimeters.  The  values  reported  in  this  work  were 
obtained  from  the  open  literature,  past  theses,  or  recent 
spectroscopic  observation  performed  here  at  AFIT.  These 
values  satisfy  the  Dunham  Equation  (Equation  I),  which  is  a 
power  series  used  to  model  the  energy  of  a  vibrating  rotator 
which  serves  as  a  model  for  the  diatomic  molecule. 

T  -  £1^  (v+l/lf  JJ  (J+l)j  (1) 

Here,  T  is  the  total  energy,  the  Yjj  's  are  known  as  the 
Dunham  coefficients,  v  is  the  vibrational  quantum  number,  and 
J  is  the  rotational  quantum  number  (15:109).  When  the 
expanded  power  series  is  truncated  so  that  j  equals  only  0  or 
1,  the  Yj.  ’s  become  the  molecular  constants  uj^  ,  u>eXe  »  Be ,  a  , 


De  ,  etc.  used  in  the  calculation  of  the  potential. 

Briefly,  the  steps  which  lead  to  the  calculation  of  the 
wavef unctions  are  as  follows.  First,  a  pair  of  classical 
turning  points  is  generated  for  each  vibrational  level. 

These  points  define  the  classical  boundary  of  the  potential 
well  and  also  "pin  down"  the  shape  of  the  well  for  future 
operations  such  as  plotting  the  potentials.  Next,  with  U(x) 
approximated  by  the  turning  points,  the  Wave  Equation  is 
solved  at  each  vibrational  level  with  J«0.  This  yields  the 
vibrational  wavefunction  which  is  used  in  connection  with 
another  vibrational  wavefunction  to  calculate  the 
Franck-Condon  Factor. 

The  calculation  of  the  potential  well  begins  with 
equation  1  which  is  rewritten  with  the  vibrational  and 
rotational  energies  split  into  seperate  terms,  Gv  and 
By,J(J+l)  respectively.  Both  terms  are  written  for  a 
vibrational  level,  and  T  is  equal  to  the  energy  (vibrational 
+  rotational)  for  that  v  level.  The  rewritten  equation  is 
shown  here  as  equation  2. 

T  -  Ew  +  Er  -  Gu  +  By J( J+l )  (2) 

Gu  becomes  the  purely  vibrational  terms  (j-0)  and  can  be 
written  as  a  power  series  in  (v  +  1/2)  (10:25  July). 

Gv  -  Z  *10  (v+1/2)1  (3) 

-Us  (v+1/2)  -u>eXe(v+l/2)2+  ... 

ByJ (J+l )  represents  the  truncated  rotational  series  (here 


j=»l)  and  the  applicable  cross  terms  (20:9). 

Bu J(  J+l  )  -  YX1  (v+1/2)1  J(J+1)  (4) 

-  Be  J(  J+l )  +cie  (v+l/2)J(J+l)  +  ... 

Equations  3  and  4  represent  the  total  vibrational  and 

rotational  energy  of  a  diatomic  molecule.  According  to  Rees, 

the  number  of  molecular  constants  required  to  equal  the 

accuracy  of  much  spectroscopic  data  is  six  (3  vibrational,  2 

rotational,  1  cross  term).  Occasionally  the  data  can  be  fit 

to  a  quadratic  and  only  five  are  required  (one  vibrational  is 

dropped)  (21:99).  The  classical  turning  points  occur  when 

the  total  energy,  T,  equals  the  potential  U.  The  minimum  and 

maximum  internuclear  separation  (turning  points)  can  be 

defined  in  terms  of  two  functions  f  and  g  (22:30).  These  two 

functions  can  in  turn  be  written  as  functions  of  the 

vibrational  and  rotational  quantum  numbers,  v  and  J,  shown  in 

equations  5a  and  5b.  ,  . 

v’+i 

f(v,J(J+l))  .  _Ti _  _ d(v+l) _ 

(2u)*  j[U-T(v+*,J(J+l))]*  (5a) 

0 

(•v 

g( v, J( J+l)  »  (2y)*  9JCJ+1) lv+1  d(v+*^  (5b) 

*  J0[U-T(v+i,J(J+l))]* 

where  u  is  the  reduced  mass  of  the  molecule.  vf  in  the  upper 
integration  limit  is  chosen  so  that  U  ■  T(v’+l/2,  J(J+1))  for 
a  fixed  value  of  J(J+1). 

If  J*0  the  rotationless  energy  states  will  be  evaluated 
leaving  only  the  vibrational  term,  Gu ,  in  the  energy 


equation.  This  simplication  will  be  applied  to  the 
calculation  of  Franck-Condon  Factors  covered  later  in  this 
section.  In  the  calculation  of  a  Franck-Condon  Factor-  the 
rotational  motion  of  the  molecule  is  ignored  and  only  the 
vibrational  wavefunction  is  used.  The  vibrational 
wavefunction  is  calculated  from  a  J-0  energy  state. 

With  J-0  equation  5a  reduces  to 

f(v)  -  [“  (Gu-  Gu.  )-*  dv'  (6) 

From  equation  2,  3T/  3J(J+1)-BV;  and  equation  5b  becomes  (for 
J-0) 

f w* 

g(v)  -  Bv,  (Gv-  Gu,  )-*  dv*  (7) 

*o 

The  lower  limit  of  integration,  vQ  ,  is  defined  by  G(vQ  )— Yqq 
or  the  electronic  energy  (22:30-31).  With  the  aid  of  a 
computer,  equations  6  and  7  can  be  solved  numerically.  Once 
f  and  g  are  obtained,  the  turning  points  can  be  calculated 
(10:25  July).  These  are  shown  in  equations  7  and  8. 

Rmin  -  |  J  +  f2]  *  -  f  (8) 

Rmax  -  [  X  +  f2  ]*  +  f  (9) 

Equations  6  and  7  also  have  physical  significance  on  the 
potential  energy  curve.  The  f  parameter  is  the  half  width  of 
the  potential  at  the  energy  level  considered  and  2g  is  the 
midpoint  of  that  energy  level.  Both  of  these  parameters  are 


required  to  accurately  position  a  vibrational  level  above  the 
internuclear  separation  axis. 

Once  the  turning  points  are  calculated,  the  first-cut 
approximation  of  the  potential  is  complete.  However,  because 
an  infinite  number  of  molecular  constants  can  not  be 
calculated  equations  3  and  4  and  therefore  the  energy 
equation  itself  must  be  truncated.  This  reduces  the  accuracy 
of  the  approximation  as  a  whole.  It  is  generally  the  upper 
states,  which  are  closer  together,  that  are  affected  the  most 
by  the  truncation.  One  method  used  to  increase  the  accuracy 
of  the  RKR  approximation  is  the  IPA  approach.  Using  the 
values  calculated  by  the  RKR  method  as  starting  values,  the 
potential  can  be  recalculated  through  an  iterative  technique 
until  it  has  the  desired  accuracy.  Yet  even  this  approach 
shows  anomalous  behavior  at  higher  energy  levels,  especially 
on  the  inner  branch  of  the  potential  energy  curve  (20:84). 

As  mentioned  before,  the  IPA  is  an  iterative  method. 

The  sequence  of  steps  which  are  iterated  are  as  follows 
(10:22  Aug): 

1)  Assume  a  function  for  the  potential  (the  RKR  value, 

for  example) 

2)  Do  the  perturbation  analysis  and  calculate  the  energy 

E 

3)  Compare  E(calculated)  with  E(measured) 

4)  If  the  two  values  differ  by  more  than  a  specified 

amount,  calculate  A U  due  to  this  difference 

5)  Use  A U  as  a  perturbation  to  the  initial  U  and 

recalulate  E( calculated ) 

Repeat  steps  one  through  five  until  the  calculated  energy 
value  is  within  a  specified  amount  of  the  measured  energy 
value.  At  each  iteration  the  potential  is  adjusted  and 


therefore  it  becomes  more  accurate  as  the  energy  values 
become  more  accurate. 

The  IPA  begins  with  the  assumption  that  a  diatomic 
molecule  can  be  modelled  as  a  rotating  vibrator.  Evidence 
that  this  is  true  can  be  found  in  the  fine  structure  of  the 
rotation  bands  (15:106).  The  rotating  vibrator  can  be 
described  by  the  radial  form  of  the  Schrodinger  Wave  Equation 
shown  in  equation  10  (27:48). 


(H0  +  Hrot  )  ¥wj(r)  -  Eyj  4*uj  (r) 


(10) 


Here,  EWj  is  the  energy  eiganvalue  (in  units  of  cm-1) 
specified  by  the  parameters  v  and  j ;  H0  is  the  Hamiltonian  of 
the  non-rotating  molecule  expanded  below  in  equation  11  (for 
units  of  cm-1). 


Hn  -  -Ti  Id* 
(4iryc ]  dr2 


UQ(r) 


(ID 


U  is  the  reduced  mass  of  the  molecule 

r  is  the  internuclear  distance 

U0(r)  is  the  potential  of  the  rotationless 

molecule 

Hrot  is  the  Hamiltonian  for  a  rotating 
molecule . 

The  exact  form  of  the  Hamiltonian  varies  for  each  coupling 
case,  but  for  a  singlet  sigma  state  it  is  (assuming  the  rigid 
rotor ) : 


rot 


_h _  fJ(J+l) 

,4irucJ  r 2 


(12) 


This  form  of  the  Schrodinger  Wave  Equation  is  then 
numerically  solved.  Initially,  the  RKR  potential  is  used  for 


Uo(r),  but  in  succeding  iterations  the  new  "perturbed"  value 
of  the  potential  is  used. 

U0(r)  -  U(00)(r)  +  AU0(r)  (13) 

■  "old  value"  +  "calculated  perturbation" 

The  solution  of  the  Schrodinger  Wave  Equation  yields  the 
energy  eigenvalue  for  the  v  and  j  being  considered.  This  is 
the  value  which  is  then  compared  to  the  measured  value  and  A  Evj 
a  Euj( measured)  -  ^(calculated)  is  calculated.  This  value  is 
then  checked  against  the  convergence  criteria.  At  this  point 
the  perturbation  analysis  is  applied.  The  new  potential 
difference  is  calculated  from  this  AE.  Normally,  AUq(t)  is 
known  and  AEvj is  the  desired  quantity. 

iE.j  -  <’:ji‘V'>in>  (u) 

The  value  on  the  right  hand  side  is  the  expectation  value. 

The  "inversion"  occurs  here  because  AEvj  is  known  and  UQ(r) 
is  the  desired  quantity,  which  must  be  chosen  to  continue  the 
evaluation.  If  the  method  is  to  converge  rapidly,  a 
well-behaved  form  of  U0(r)  must  be  used.  Good  results,  with 
the  trade  off  of  longer  compiling  times,  have  been  acheived 
with  a  form  of  AUQ(r)  which  is  a  summation  of  a  Legendre 
Polynomial  damped  by  an  exponential  function  (27:50).  It  is 
this  form  of  the  potential  which  is  used  in  the  programs 
presented  in  this  paper. 

With  the  calculation  of  AUg(r)  the  Hamiltonian  becomes. 
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H  -  H0  +  Hrot  +  AUq (r)  (15) 

This  Hamiltonian  is  now  used  in  the  Schrodinger  Wave  Equation 
which  is  solved  numerically  to  get  a  new  value  of  Evj 
(calculated),  and  the  iteration  procedes  as  usual.  When  the 
convergence  criteria  is  met,  the  process  is  stopped.  With 
the  improved  potential,  more  accurate  wavef unctions  can  be 
calculated,  and  from  them  more  accurate  Franck-Condon 
Factors . 

Franck-Condon  Factors 

The  desired  outcome  from  these  calculations  is  a  set  of 
numbers  between  zero  and  one,  one  for  each  possible 
transition  between  vibrational  levels  of  the  two  states.  The 
closer  the  number  is  to  one  the  larger  the  probability  of 
that  transition  occurring,  and  therefore  the  stronger  its 
intensity  will  be.  Because  of  the  normalization  of  the 
wavef unctions ,  the  total  of  all  Franck-Condon  Factors  for 
transitions  from  or  to  a  specific  vibrational  level  should  be 
one.  The  process  of  calculating  Franck-Condon  Factors  comes 
full  circle  when  a  spectroscopist  uses  a  set  of  "raw"  factors 
to  justify  the  labelling  of  transitions  in  a  spectrum  used  to 
calculate  the  constants  which  in  turn  are  used  to  calculate 
the  potentials  that  the  factors  are  derived  from. 

The  theory  behind  the  math  was  originally  proposed  by 
James  Franck  in  1925  to  explain  the  dissociation  of  iodine 
vapor  (13:21).  Edward  Condon  expanded  Franck’s  principle  and 
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applied  it  to  intensity  distributions  in  1926.  In  1928, 
Condon  updated  the  principle  to  incorporate  principles  from 
the  then  new  quantum  physics  (6:32).  It  is  in  this  form  that 
the  principle  now  stands.  The  underlying  concept  of  the 
Franck-Condon  Principle  is  really  quite  simple.  Because  the 
mass  difference  between  the  nuclear  material  and  the 
electrons  is  quite  large,  when  the  electron  jumps  from  one 
electron  state  to  another  the  transition  takes  place  "so 
rapidly  in  comparison  to  the  vibrational  motion  that 
immediately  afterwards  the  nuclei  still  have  very  nearly  the 
same  relative  position  and  velocity  as  before  the  'jump'" 
(15:194).  Or  stated  another  way,  the  transition  will  be 
vertically  upward  or  downward  onto  the  potential  energy  curve 
of  the  second  state.  The  strength  of  the  transition  will 
then  be  determined  by  how  much  time  the  electron  spends  at 
both  the  before  and  after  location.  In  other  words,  the 
stronger  transitions  occur  between  locations  of  maximum 
probability  for  the  electron.  Classically,  the  electron 
spends  most  of  its  time  at  the  minimum  and  maximum 
intersections  of  the  vibrational  level  with  the  potential,  as 
seen  in  Figure  I.  The  most  probable  (most  intense) 
transitions  would  occur  where  the  upper  turning  point  is 
vertically  above  the  lower  point  as  in  transitions  A  and  B. 
But  in  transition  C,  the  upper  turning  point  is  vertically 
above  the  midpoint  of  the  vibrational  level  of  the  lower 
state.  Therefore  C  will  be  a  weak  (less  probable) 
transition.  To  determine  strengths  of  transitions  quantum 


mechanically,  the  wavef unctions  for  each  vibrational  level 
must  be  studied.  As  in  the  classical  case,  the  most  probable 
transition  will  occur  when  a  maximum  of  the  upper  state  lies 
vertically  above  a  maximum  of  the  lower  state.  As  higher 
vibrational  levels  are  considered,  the  points  of  maximum 
probability  for  the  position  of  the  electron  shift  towards 
the  classical  turning  points.  But  because  quantum  mechanics 
is  more  concerned  with  the  probability  of  electronic  position 
rather  than  its  actual  location  the  value  of  the 
wavef unctions  must  be  integrated  over  the  entire  vibrational 
level.  The  amount  of  overlap  of  the  two  levels  of  a 
transition  is  what  then  becomes  important.  There  has  even 
been  some  discussion  in  the  literature  (23)  that  mention  of 
"vertical"  transitions  should  be  dropped  from  any  description 
of  the  Franck-Condon  Principle  on  a  quantum  mechanical  level. 
In  other  words,  the  only  statement  necessary  is  one  of  the 
overlap  integral.  From  this  point  of  view,  the  transitions 
are  no  longer  required  to  be  vertical.  The  manner  of  overlap 
is  also  important.  For  maximum  probability,  the  maximum  of 
the  upper  and  lower  states  should  occur  at  the  same 
separation  of  the  two  molecules.  A  minimum  on  one  level  will 
cancel  out  a  maximum  on  the  other.  The  product  of  the  two 
vibrational  wavef unctions  integrated  across  their  overlap  and 
squared  is  the  Franck-Condon  Factor  (18). 

The  Franck-Condon  Factor  can  also  be  derived  from  a 


quantum  mechanical  statement  of  the  transition  moment  between 
an  upper  state  v'  and  a  lower  state  v".  The  probability  of 


the  transition  occurring  (known  as  the  Einstein  Coefficient) 
is  directly  proportional  to  the  square  of  the  transition 
moment  (also  known  as  the  electric  moment).  The  transition 
moment  is  a  function  of  the  dipole  moment  (M);  a  vector  whose 
components  are  the  sum  of  each  charge  multiplied  by  its 
distance  from  the  origin  (15:199).  The  transition  moment  is 
shown  in  equation  16. 

R  =  f*? ' #M4'"  dx  , 


By  referring  to  the  Born-Oppenheimer  Approximations  and 
ignoring  rotational  effects,  the  two  molecular  wavef unctions 
can  be  split  into  a  vibrational  and  an  electronic  part.  That 
is ,  4*  =  4^  *  4^;  this  is  possible  because  the  nuclear  motion  is 
much  slower  than  electron  motion  in  the  molecule  and  thus  the 
wavefunction  describing  a  particular  vibrational  level  of  an 
electronic  state  can  be  broken  down  into  the  product  of  two 
wavef unctions ,  one  purely  for  the  electronic  state,  the  other 
for  the  vibrational  level  under  consideration.  Each  of  the 
separated  wavef unctions  is  now  dependent  upon  independent 
coordinate  systems,  one  for  the  electrons,  the  other  for  the 
nuclear  components  (17:231).  Similiarly,  the  dipole  moment, 
M,  can  be  broken  down  into  the  sum  of  a  nuclear  part  based  on 
nuclear  coordinates  and  an  electronic  part  based  on  the 
electron  coordinate  system,  M  ■  Me  +  Mn  (15).  Both  expanded 
equations  can  then  be  substituted  back  in  equation  16. 

Because  the  vibrational  wavefunction  is  real,  its  complex 
conjugate  is  identical  to  the  original  wavefunction.  This 


fact  and  the  expanded  form  of  16  allow  the  transition  moment 
to  also  be  written  as  a  sum  of  an  electronic  and  nuclear  part 
as  in  equation  17. 
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Looking  at  the  second  term  and  remembering  that  the  two 
coordinates,  electronic  and  nuclear,  are  independent  of  one 
another,  the  integral  can  be  rewritten  as: 

|m  rr'  dx  frV'  dx  as) 

J  n  V  V  n  J  e  e  e  v  ' 

The  second  integral  in  equation  18  has  a  value  of  either  zero 

or  one.  The  set  of  electronic  wavef unctions  is  orthonormal 

and  the  integration  of  their  product  from  plus  to  minus 

infinity,  will  have  a  value  of  one  if  VJ  and  zero  if  'fe'  * 

4*"  If  the  former  were  true  no  transition  would  have 
e 

occurred,  therefore  the  value  of  the  integral  and  further, 
the  entire  second  term  is  zero. 

The  first  term  can  also  be  split  into  the  product  of  a 
nuclear  and  electronic  part.  The  equation  for  the  transition 
moment  reduces  to  the  form  in  equation  19. 

R  =  [  4* '  4* "  dx  [m  4"*4'"  dx  (19) 

J  v  v  n  J  e  e  e  e 

The  second  integral  is  known  as  the  electronic  transition 
moment,  Re .  This  factor  does  depend  on  internuclear  distance 


but  here  the  Franck-Condon  Principle  is  brought  into  play.  A 
restatement  of  the  principle  is  that  Re  varies  slowly  with 
internuclear  distance  and  therefore  may  be  replaced  by  a 


constant  value  TJe.  The  constant  value  can  then  be  pulled  out 
of  the  Integral  and  the  transition  moment  becomes: 

R  *  R  Iff"  dT  (20) 

ej  v  v  n 

The  remaining  integral  should  be  recognized  as  the  overlap 
integral,  and  its  square  as  the  Franck-Condon  Factor  (10:26 
Sept).  The  approximation  made  earlier,  of  the  Franck-Condon 
Factor  being  proportional  to  the  transition  probability  is 
valid  only  as  long  as  Rg ,  which  is  also  known  as  the  R 
centroid,  remains  constant.  There  are  some  who  believe  this 
approximation  is  not  accurate  enough  and  attempts  have  been 
made  to  calculate  directly  the  fluctuation  of  Rg(r)  with  r 
(4:77,  18:95).  It  was  hoped  that  gains  would  be  made  in 
three  areas:  a  form  of  Re(r)  vs  r  would  be  found, 
intensities  of  unmeasured  bands  could  be  calculated,  and 
experimental  data  obtained  from  several  measurements  could  bs 
more  easily  smoothed  (16:20).  This  area  of  study  known  as 
the  R  Centroid  Approach  has  generally  met  with  good  but 
limited  results  (the  method  does  not  accurately  describe 
highly  quantum  transitions,  i.e.  transitions  not  classically 
allowed).  This  approach  has  been  criticized  by  Thomas  James 
(16)  who  concluded  that  the  R  Centroid  Approach  may  be  "an 
unsatisfactory  method  of  treating  experimental  data  in  many 
cases" . 

Returning  to  the  overlap  integral,  because  the 
vibrational  wavef unctions  represent  different  electronic 
states  they  are  not  an  orthonormal  set,  and  the  integral  will 


have  a  value  other  than  zero  or  one.  Each  wavefunction  is 
normalized  however;  and  has  a  maximum  value  of  one. 

Therefore  the  maximum  value  for  the  integral  is  one  and  in 
general,  it  will  vary  between  zero  and  one.  The  evaluation 
of  this  integral  is  the  key  to  a  computer  calculation  of 
Franck-Condon  Factors.  A  compromise  must  be  reached  between 
accuracy  of  the  method  and  the  computer  time  required  to 
implement  it.  A  Simpson's  Rule  approximation  to  the  overlap 
integral  was  sufficient  for  the  calculations  in  this  work 
(see  equation  21). 

{Vi  dX  **  (2i) 

+(4f1+2.Atit3)] 

Here  the  Al^'s  represent  the  product  of  the  upper  and  lower 
wavefunction  values  at  location  *if.  In  other  words,  ATj.Af'i 
*  A4*i  .  Simpson's  Rule  improved  accuracy  over  the  previous 
method,  a  histogram  integral  approximation,  but  did  not 
noticably  increase  the  run  time  of  the  program. 


III.  COMPUTER  PROGRAMS 
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Program  RKR-IPA 

The  potential  modelling  program,  RKR-IPA,  was  written  by 
Dr  C.R.  Vidal  in  Fortran  IV  to  run  on  a  Cray-1  computer.  The 
program  was  brought  to  AFIT  and  altered  to  run  on  available 
systems  in  1983  by  Capt  Lyle  Rutger  as  a  part  of  his  thesis 
project.  The  program  was  changed  to  Fortran  V  and  the  arrays 
were  decreased  in  size  to  meet  memory  requirements  for  the 
CDC  Cyber  located  at  Wright-Patter son  AFB.  Other 
modifications  were  made  to  the  program  to  update  the  code  and 
conform  its  output  to  Rutger's  uses.  The  program  calculates 
an  RKR  potential  and  contains  an  optional  IPA  routine  which 
increases  the  accuracy  of  the  potential  approximation.  Input 
to  the  program  is  in  the  form  of  a  data  file  containing  all 
the  necessary  information  to  run  the  program.  The  RKR-IPA 
program  is  not  currently  set  up  to  run  interactively.  Simple 
modifications  would  allow  this  to  be  done,  but  because  of  the 
large  amount  of  output  generated  by  the  program  it  has  been 
kept  in  its  batch  configuration.  The  RKR-IPA  program  also 
produces  two  computer  files.  The  first  is  a  file  referred  to 
as  "POINTS"  in  the  program.  This  file  contains  the  turning 
points  for  each  vibrational  level  of  the  potential.  POINTS 
is  used  as  input  to  the  plotting  program  (CRVPLT)  written  by 
Pow.  The  second  file  created  is  called  "WAVE"  and  is  of 
primary  importance  to  the  calculation  of  the  Franck-Condon 
Factors.  WAVE  can  be  broken  down  into  several  sections,  its 
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first  two  lines  cover  general  information  about  the  potential 
as  a  whole.  Line  one  shows  the  highest  vibrational  level  for 
which  a  wavefunction  has  been  calculated;  line  two  contains 
the  minimum  and  maximum  internuclear  distances  over  which  the 
potential  has  been  calculated.  The  third  number  on  line  two 
is  the  number  of  points  at  which  the  wavefunction  has  been 
evaluated.  The  remainder  of  the  file  is  composed  of 
tabulated  values  of  the  wavef unctions  for  each  vibrational 
level  in  the  following  format.  The  first  line  of  each 
section  contains  information  about  a  specific  level  and  has 
three  numbers  on  it.  The  first  is  the  vibrational  level  of 
the  section  (starting  at  zero),  the  second  is  the  array 
location  of  the  first  non-zero  (greater  than  1  E-8)  value  of 
the  wavefunction,  and  the  third  is  the  array  location  of  the 
last  non-zero  value  of  the  wavefunction.  These  last  two 
numbers  are  referred  to  as  the  start/stop  indices  in  the 
program  documentation.  The  next  N  lines,  where  N  is  the 
number  of  points  given  in  line  two,  contains  the  tabulated 
values  of  the  wavefunction.  One  value  for  each  array 
location  which  corresponds  to  an  internuclear  distance 
through  a  scaling  factor.  This  same  pattern,  one  line  of 
information  followed  by  N  lines  of  wavefunction  values,  is 
repeated  for  the  rest  of  the  file.  The  final  section  is  for 
the  highest  vibrational  level  referred  to  in  line  one  of  the 
f  ile . 

An  example  of  the  input  data  file  is  shown  in  Figure  II. 
It  contains  all  of  the  information  needed  by  the  program  to 
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model  the  potential  and  allows  the  selection  of  all  available 


options.  It  also  indicates  how  many  vibrational  levels  will 
be  calculated.  The  information  is  presented  in  the  following 
format  (the  numbers  correspond  to  the  lines  indicated  in  the 
figure) : 

1.  The  first  line  should  have  a  one  in  column  one  to 
insure  the  program  will  execute  rather  than  just  list.  Using 
the  /NOSEQ  command  in  the  jobfile  used  to  submit  the  batch 
job  will  insure  that  the  fijrst  column  of  the  data  file  is 
processed  by  the  computer.  Without  this  command,  the 
computer  will  assume  the  first  several  spaces  contain  line 
numbers  and  it  will  not  process  any  information  contained  in 
these  spaces. 

2.  The  second  and  third  lines  indicate  the  molecule  and 
state  to  be  investigated. 

3.  The  fourth  line  contains  the  mass  scale  to  be  used 
and  the  masses  of  the  two  atoms.  A  one  in  column  four 
indicates  the  masses  are  in  the  Carbon-12  scale,  a  two 
indicates  oxygen-16.  The  masses  are  entered  next  with  a 
2D16.9  format.  The  reduced  mass  of  the  molecule  may  be  used 
instead  of  the  individual  masses  by  entering  it  with  the  same 
format  in  place  of  the  first  mass  and  leaving  the  second  mass 
blank . 

4.  The  number  of  Dunham  Coefficients  presented  is 
covered  in  the  next  line.  Format  is  415.  The  first  number 


is  the  number  of  lines  containing  coefficients,  the  remaining 
numbers  give  the  number  of  coefficients  in  each  line.  For 


example,  a  common  sequence  is  3,  3,  2,  1.  This  would  be 
interpreted  as  three  lines  of  coefficients  follow,  three 
coefficients  on  the  first  line  (Ta  ,we,  <oeXe),  two  on  the 
second  (Be,ae),  and  one  on  the  third  (De). 

5.  Enter  the  constants  as  described  in  point  4.  Format 
is  D18.9. 

6.  The  next  line  after  the  final  Dunham  Coef ficient(s) 
contains  the  discociation  energy  and  the  term  energy.  Format 
is  2D15.8. 

7.  Next,  indicate  the  highest  vibrational  level  to  i,e 
investigated,  and  follow  by  the  incremental  step  size  used  in 
calculating  the  turning  points.  Format  is  2F6.2. 

8.  If  the  values  of  the  Klein  Action  Integral  are  to  be 
printed,  the  first  column  of  the  next  line  should  contain  a 
one.  A  zero  will  skip  printing  the  values. 

9.  If  Vidal's  extension  routine  is  to  be  used  to  extend 
the  curve  past  the  normal  range  of  the  program,  an  extended 
internuclear  range  must  be  defined.  This  range  should  be 
large  enough  to  contain  the  minimum  and  maximum  internuclear 
separation  represented  by  the  extended  potential  well.  The 
format  for  these  two  separations  is  2F6.2.  The  final  number 
on  the  line  is  the  number  of  extra  data  points  used  for  the 
extension.  This  number  must  be  between  0  and  25,  with  a 
format  of  12. 

10.  If  the  last  number  on  the  previous  line  is  anything 
other  than  zero,  the  data  points  are  entered  on  this  line  as 
described  in  Rutger's  thesis  (22:B-9).  Otherwise  the  line  is 


left  blank. 


11.  The  next  line  should  contain  two  zeroes,  format 
212.  This  line  is  a  hold  over  from  when  the  program  punched 
an  output  deck.  Since  Rutger  changed  the  PUNCH  commands  to 
comment  lines  ,  that  option  is  no  longer  available  without 
modifying  the  program.  A  zero  suppresses  punching  the  card 
decks . 

12.  If  the  IPA  routine  is  desired  this  line  should 
contain  a  1  in  column  one.  A  zero  will  skip  the  IPA  and 
repeat  the  RKR  analysis  on  the  next  set  of  data.  If  no 
further  data  is  given,  this  line  should  be  followed  by  a  zero 
on  a  seperate  line. 

Note:  The  following  data  lines  are  needed  only  if  the  IPA 

option  is  selected. 

13.  A  one  in  column  four  on  this  line  will  print  every 
iteration  of  the  function  SCHR  (the  function  subroutine  which 
solves  the  Schrodinger  Wave  Equation),  a  zero  will  suppress 
printing.  A  one  in  column  eight  will  print  the  wavefunction 
value  at  every  Nth  point,  where  N  is  given  by  the  third 
number  on  the  line.  A  zero  in  column  eight  will  not  print 
the  values.  Columns  nine  through  twelve  are  for  N  which 
defines  how  many  values  of  the  wavefunction  will  be  printed. 
The  next  four  columns  are  for  the  number  of  iterations  SCHR 
will  complete  in  an  attempt  to  find  a  converging  solution. 
Finally,  enter  the  convergence  criteria  for  the  Schrodinger 
Wave  Equation  with  a  D10.0  format. 

14.  Next,  define  the  range  of  internuclear  separations 
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over  which  the  Schrodinger  equation  will  be  analyzed.  Format 
is  2F10.1,  the  values  are  in  units  of  angstroms.  The  first 
number  represents  the  first  position  of  the  wavefunction 
array  (minimum  internuclear  separation)  and  the  second  number 
represents  the  last  position  of  the  array  (maximum 
separation).  Beginning  in  column  21  enter  the  number  of 
elements  in  the  wavefunction  arrays,  (that  is,  the  number  of 
points  at  which  the  wavefunction  will  be  evaluated)  with  an 
15  format.  The  maximum  number  is  2403.  Beginning  in  column 
26  with  an  15  format  enter  the  number  of  iterations  of  the 
IPA  routine  to  be  performed.  Starting  in  column  31  enter  the 
number  of  terms  to  be  used  in  the  Legendre  polynomial,  format 
15.  This  is  part  of  the  chosen  function  used  for  AUo(r). 
Next,  place  a  zero  in  column  36.  This  is  another  hold  over 
from  output  punched  on  card  decks.  Finally,  in  column  40 
should  come  a  zero  if  only  the  first  and  last  iterations  of 
IPA  are  to  be  printed,  or  a  one  if  all  iterations  should  be 
printed . 

15.  The  next  two  lines  contain  eleven  sets  of  three 
number  groups  which  define  the  range  of  vibrational  and 
rotational  levels  to  be  studied.  The  first  line  contains  six 
sets,  the  second  five.  Each  group  has  a  format  for  the  three 
numbers  of:  12,  13,  F7.0.  The  order  within  the  group  is, 
minimum  V+l  level,  maximum  V+l  level,  J(J+1)  value.  Note 
that  the  range  of  vibrational  levels  studied  for  a  particular 
J(J+1)  value  are  entered  as  V+l  values,  not  V.  The  last 
group  of  three  numbers  on  these  two  lines  must  have  J(J+1)  - 
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0  so  that  the  rotationless  potentials  will  be  evaluated.  If 
the  number  of  vibrational  levels  calculated  and  read  to  the 


file  WAVE  is  to  be  changed,  the  new  V  level  must  be  shown 
both  in  the  final  group  of  numbers  in  these  two  lines  and  in 
the  first  number  described  in  point  7.  The  maximum  v+1  level 
for  the  rotationless  group  (J(J+1)  should  be  the  highest 
vibrational  level  called  for  in  these  two  lines.  This  will 
ensure  that  when  the  file  WAVE  is  created  it  will  contain 
exactly  the  number  of  vibrational  levels  desired,  and  no 
extra  ones. 

16.  The  last  line  must  contain  a  zero  to  halt  execution 
of  the  program. 

The  most  frequent  changes  to  the  data  file  for  a 
particular  electronic  state  are  to  the  number  of  v  levels 
calculated  and  the  number  of  points  at  which  the 
wavef unctions  are  evaluated.  The  lines  which  must  be  altered 
are  covered  under  points  14  and  7  for  the  former  and  point  13 
for  the  latter.  For  a  thorough  explanation  of  the  variables 
represented  in  the  data  file  and  the  calculations  done  in  the 
RKR-IPA  program  please  refer  to  Rutger's  thesis,  Appendix  B 
(22:B-1).  For  a  program  listing  of  the  RKR-IPA  code  refer  to 
Pow's  thesis,  Appendix  B. 

Program  FCFACT 

Program  FCFACT  was  written  by  Capt  Pow  in  Fortran  V  to 
calculate  Franck-Condon  Factors  from  wavefunction  values 
calculated  by  the  RKR-IPA  program.  FCFACT  uses  the  files 
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created  by  two  RKR-IPA  runs  as  input,  one  for  the  upper 
electronic  state  and  the  second  for  the  lower.  Before  any 
calculations  can  take  place,  the  potentials  must  be  checked 
to  insure  proper  overlap  and  scaling.  The  scaling  factor 
(angstroms  per  array  location)  is  calculated  in  line  720  for 
the  upper  molecular  state  and  in  line  730  for  the  lower.  If 
the  scaling  of  the  two  states  differs  by  more  than  a 
specified  amount,  the  lower  state  is  rescaled  to  match  the 
upper  state.  The  method  of  rescaling  is  described  in  Pow's 
thesis  (20:44-47).  Once  the  scaling  has  been  adjusted,  the 
overlap  of  the  two  states  is  checked,  and  corrected  if 
necessary.  Once  again,  the  values  of  the  lower  wavefunction 
are  shifted  to  match  the  upper.  The  purpose  of  this 
operation  is  to  insure  that  matching  array  locations 
represent  the  same  internuclear  separation.  The  shift 
subroutine  is  described  on  page  47  of  Pow's  thesis  (20:47). 
Finally,  the  wavef unctions  of  both  states  are  normalized  so 
that  the  sum  of  the  squares  of  the  wavefunction  values  for 
one  vibrational  level  equals  one.  Pow  discusses  the 
normalization  subroutine  on  page  48  of  his  thesis  (20:48). 
After  these  three  processes  have  been  completed,  the 
wavef unctions  are  properly  prepared  to  begin  calculation  of 
the  Franck-Condon  Factors. 

The  calculation  of  Franck-Condon  Factors  involves  a 
multiplication  of  the  values  of  two  wavef unctions ,  one  from 
the  upper  state  of  the  transition  and  one  from  the  lower 
state,  plus  an  integration  across  the  region  where  the  two 


states  overlap.  The  result  is  then  squared  to  give  the 
desired  quantity.  This  process  is  carried  out  for  each 
possible  combination  of  upper  and  lower  vibrational  level. 
Originally  the  integral  was  approximated  by  a  histogram.  The 
values  at  each  array  location  of  the  two  wavef unctions  were 
multiplied  together  and  the  individual  products  were  summed. 
The  start  and  stop  locations  of  the  region  of  overlap  were 
determined  from  the  start/stop  indices  before  beginning  the 
integration.  Once  the  final  pair  of  values  was  multiplied 
and  summed,  the  number  was  squared  and  stored  in  an  array 
called  FACTOR.  Once  all  of  the  Fr anck-Condon  Factors  were 
calculated  the  array  was  printed  out.  FCFACT  has  been 
revised  for  this  thesis,  to  use  a  Simpson's  Rule 
approximation  of  the  integral.  The  loop  which  integrates  and 
calculates  the  Franck-Condon  Factors  falls  between  lines  870 
and  1100  in  the  program.  First,  the  tabulated  values  of  the 
wavef unctions  for  the  two  vibrational  levels  of  interest  are 
multiplied  at  matching  array  locations  and  the  product  stored 
in  an  array  named  WFXN.  The  Simpson's  Rule  approximation  is 
carried  out  on  this  array  and  makes  use  of  a  conventional 
code  (24).  Because  of  the  speed  of  the  approximation,  the 
integration  is  carried  out  over  the  entire  product  array 
which  has  the  same  length  as  the  shorter  of  the  two 
wavefunction  arrays,  rather  than  just  over  the  non-zero 
overlap  of  the  two  wavef unctions .  Because  of  the  tabulated 
nature  of  the  wavefunction  values,  the  Ax  required  in  the 
Simpson's  Rule  approximation  (see  equation  21)  is  not  the 
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scaling  factor  SCALUP,  but  one.  This  is  because  the 
approximation  steps  through  the  array  using  each  element  in 
the  calculation.  As  before,  the  final  value  of  the  integral 
is  squared  and  stored  in  the  FACTOR  array,  to  be  printed  out 
after  all  values  have  been  calculated.  The  Simpson's  Rule 
approximation  was  verified  by  comparison  with  a  set  of  values 
known  to  be  correct.  These  values  were  calculated  by  the 
earlier  histogram  version  of  FCFACT  from  a  set  of  test  data 
initially  used  to  insure  the  program  was  producing  accurate 
numbers . 

Several  changes  must  be  made  to  program  FCFACT  in  order 
to  prepare  it  to  calculate  the  Franck-Condon  Factors  for  a 
specific  molecular  transition.  The  arrays  in  the  program 
must  be  dimensioned  to  hold  the  required  number  of  elements. 
The  first  dimension  is  for  the  vibrational  levels  and  begins 
with  zero,  the  second  represents  the  number  of  points  and 
begins  with  one.  The  simplest  way  to  avoid  having  to 
re-dimension  all  the  arrays  before  every  run  is  to  find  a 
reasonable  combination  of  V  levels  and  points  which  is  larger 
than  or  equal  to  the  values  in  the  runs  to  be  made.  These 
numbers  need  not  be  changed  until  a  larger  run  is  required. 
All  of  the  changes  described  below  must  be  made  internal  to 
the  program  and  can  be  done  using  the  line  or  screen  editor 
available  on  the  operating  system  of  the  CDC  Cyber.  The  two 
wavefunction  arrays  must  be  identically  dimensioned  for  the 
program  to  run.  Once  the  arrays  are  correctly  dimensioned, 
the  corresponding  DATA  statements  must  be  updated  to  insure 
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that  each  array  location  will  be  initialized  to  zero.  The 
large  numbers  in  the  DATA  statements  represent  the  total 
number  of  elements  in  the  array,  calculated  by  multiplying 
the  first  array  dimension  (V+l)  by  the  second  (number  of 
points).  These  DIMENSION  and  DATA  changes  must  be  made  the 
entire  length  of  the  code,  wherever  the  appropriate 
statements  or  numbers  are  found,  the  subroutines  also  contain 
several  arrays  which  must  be  changed.  Next,  the  names  of  the 
files  opened  in  lines  380  and  390  must  be  changed  to  the 
names  of  the  files  containing  the  wavefunction  values  for  the 
two  states  in  the  transition  of  interest.  The  lower  state 
must  always  be  opened  under  unit  7  and  the  upper  under  unit 
8.  Finally,  the  header  for  the  matrix  of  Franck-Condon 
Factors  must  be  changed  to  show  the  correct  transition  and 
any  other  pertinant  information.  Program  FCFACT  in  its 
present  form  can  be  run  interactively,  if  the  local  directory 
containing  the  compiled  version  of  the  program  also  contains 
the  two  input  data  files.  An  interactive  run  can  be  a  useful 
debugging  tool  if  available  options  on  the  Fortran  V  compiler 
are  used.  Compiling  the  program  with  these  debugging  options 
choosen  and  then  running  it  interactively  can  point  out 
errors  more  quickly  than  waiting  for  the  paper  output  from  a 
batch  job.  The  code  for  the  updated  portion  of  FCFACT  is 
contained  in  Appendix  A  of  this  thesis. 


IV  Results  and  Analysis 

Franck-Condon  Factors  were  calculated  for  six  molecules, 
lead  oxide  and  five  halogen  or  interhalogen  compounds  (Br2 , 
BrCl,  BrF,  IC1 ,  IF).  The  molecular  constants  for  lead  oxide 
were  calculated  from  data  collected  during  spectroscopic 
observation  in  the  lab  at  AFIT.  The  potentials  and 
wavefunction  values  were  calculated  for  nine  states  of  PbO 
(a,  b.  A,  B,  C,  C’,  D,  E,  X)  and  Franck-Condon  Factors  were 
calculated  for  the  transition  from  each  state  to  ground,  plus 
the  D  to  a  and  A  to  a  transitions.  The  B  to  X  transitions  of 
this  molecule  were  used  as  the  test  data  for  initial 
verification  of  the  programs  and  all  modifications.  The 
molecular  constants  for  the  halogens  were  drawn  from  the  open 
literature.  The  number  of  potentials,  wavef unctions ,  and 
transitions  vary  for  each  halogen.  The  results  for  each 
molecule  follow,  with  lead  oxide  first  and  then  the  halogens 
in  alphabetical  order. 

Lead  Oxide 

Lead  oxide  was  thoroughly  investigated  in  a  previous 
thesis  (20)  using  constants  calculated  by  Rutger  (22:54)  and 
from  the  open  literature.  RKR-IPA  and  FCFACT  were  run  using 
those  values  to  insure  valid  copies  of  the  programs  had  been 
obtained.  The  test  run  used  the  data  for  the  B  and  X  states 
of  PbO  and  calculated  15  vibrational  levels  for  each.  The 
results  were  verified  by  comparison  wih  the  values  Pow 


obtained.  Next,  the  programs  were  rerun  using  the  same  data, 
calculating  the  wavefunction  values  and  Pranck-Condon  Factors 
for  24  vibrational  levels.  The  programs  were  again  run  for 
24  levels  using  the  new  molecular  constants  calculated  by 
Capt  Jeff  Duray  and  presented  in  his  thesis  (11).  The 
Franck-Condon  Factors  from  the  two  sets  of  data  were  compared 
for  each  transition  with  the  largest  change  being  seen  in  the 
higher  vibrational  levels.  This  was  expected,  as  the  purpose 
of  the  new  investigation  of  lead  oxide  was  to  increase  the 
accuracy  of  the  constants  for  the  higher  levels.  The  Condon 
Parabola  is  a  characteristic  feature  of  tables  of 
Franck-Condon  Factors,  shown  by  two  lines  of  intensity  maxima 
in  the  table.  These  are  due  to  the  fact  that  there  will  be 
two  values  of  maximum  probability  for  each  level  of  the 
beginning  state  in  the  transition,  one  for  each  turning  point 
of  that  level,  see  Figure  I.  Theoretically,  these  two  lines 
will  form  a  parabola  (15:196-198).  The  Condon  Parabola  for 
the  new  data  shows  only  minor  deviations  from  the  older 
parabola  and  these  are  for  v>10.  For  v<10  the  two  parabolas 
coincide  exactly.  The  two  sets  of  Franck-Condon  Factors,  for 
the  region  of  the  parabola,  are  shown  in  Table  1.  The  first 
value  is  the  value  calculated  by  Pow  using  a  simple  Histogram 
approach.  The  second  is  the  value  calculated  using  Duray's 
molecular  constants  and  the  same  Histogram  approach.  The 
third  entry  is  the  absolute  difference  between  Franck-Condon 
Factors  calculated  using  in  one  case  the  simple  Histogram 
approach  and  in  the  second  case  the  more  accurate  Simpson's 
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Rule  approximation  to  the  integral  (both  with  Duray's  data). 
The  PbO  B  to  X  transition  was  again  used  as  the  test  case  for 
this  modification  to  the  program  FCFACT. 

All  of  the  states  for  PbO  presented  in  this  thesis  were 
calculated  from  the  more  recent  data,  which  is  presented  in 
Table  2.  The  values  of  the  wavef unctions ,  except  for  the  E 
state,  were  calculated  for  a  total  of  25  levels  (beginning  at 
zero)  and  evaluated  at  301  points  on  each  vibrational  level. 
It  is  important  to  stay  within  the  accurate  portion  of  the 
potential  energy  curve  when  extending  the  number  of  levels 
calculated.  The  potential  approximation  begins  to  break  down 
near  the  dissociation  energy  and  special  procedures  need  to 
be  undertaken  to  maintain  the  desired  level  of  accuracy  in 
this  region.  For  all  states  of  PbO  except  the  E  state,  level 
24  is  within  the  reliable  region.  The  Franck-Condon  Factors 
for  each  of  the  excited  states  to  the  ground  state  follow  in 
Tables  3  through  10.  The  factors  for  the  two  transitions 
between  excited  states  calculated,  A  to  a  and  D  to  a,  are 
included  in  Tables  11  and  12  respectively.  The  E  state  has 
the  highest  potential  energy  of  all  the  PbO  excited  states 
and  is  very  wide  and  shallow.  The  state  is  difficult  to 
calculate  and  the  maximum  vibrational  level  accurately 
attainable  is  15.  Because  of  the  overlap  between  the  E  and 
the  X  states  and  because  of  the  requirements  of  the  program 
FCFACT,  a  minimum  of  351  points  were  required  to  give  a  valid 
calculation  of  the  Franck-Condon  Factors  for  the  transitions 
between  these  two  states. 
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TABLE  3 


Franck-Condon  Factors  for  PbO  Small  A  to  X  Transitions 
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Values  in  the  table  are  (FCF  x  1000) 


2 

067 

029 

010 

003 

000 

000 

000 

000 

000 

000 

000 

000 

3 

136 

108 

060 

025 

002 

000 

000 

000 

000 

000 

000 

000 

4 

034 

108 

130 

095 

018 

005 

001 

000 

000 

000 

000 

000 

5 

038 

000 

054 

120 

078 

035 

012 

003 

001 

000 

000 

000 

.i 

6 

040 

071 

017 

010 

128 

106 

058 

023 

007 

001 

000 

000 

h  *  •  * 

c 

7 

032 

003 

060 

053 

035 

109 

125 

084 

038 

013 

003 

001 

8 

015 

056 

008 

023 

026 

004 

071 

126 

108 

059 

022 

006 

9 

046 

001 

039 

039 

052 

057 

003 

031 

107 

125 

082 

036 

10 

001 

037 

023 

009 

Oil 

020 

068 

025 

005 

075 

127 

105 

11 

027 

023 

009 

044 

037 

037 

001 

053 

052 

001 

039 

114 

12 

031 

003 

039 

001 

021 

010 

051 

007 

025 

066 

017 

012 

A 

13 

001 

034 

005 

027 

Oil 

040 

000 

041 

029 

004 

059 

041 

W 

14 

017 

018 

013 

026 

036 

000 

039 

013 

018 

047 

001 

038 

Pr,-T 

15 

031 

001 

032 

000 

001 

032 

Oil 

020 

032 

002 

047 

016 

V-\  *• 

16 

010 

022 

008 

023 

020 

016 

013 

029 

003 

041 

003 

031 

>.V." 

17 

001 

026 

004 

024 

026 

003 

031 

000 

035 

002 

032 

019 

18 

019 

005 

025 

002 

001 

028 

002 

028 

005 

025 

016 

014 

.  • .  * 

•  v 

A 

19 

026 

003 

020 

009 

013 

013 

015 

015 

013 

020 

008 

030 

20 

010 

020 

002 

025 

025 

001 

025 

002 

027 

001 

030 

000 

rz 

21 

000 

023 

005 

015 

007 

020 

004 

023 

002 

025 

003 

026 

.  ■/. .  • 

22 

008 

008 

020 

000 

001 

020 

006 

016 

010 

014 

012 

015 

23 

021 

000 

020 

007 

017 

002 

022 

000 

023 

001 

024 

002 

24 

021 

008 

006 

020 

020 

005 

013 

012 

007 

018 

003 

022 

V 

TABLE  4 


Fr anck-Condon  Factors  for  PbO  Small  B  to  X  Transitions 


2  3  4  5  6  7 


10  11 


0  016  066  138  192  201  167  113  063  029  011  003  001 


1  073  164  151  056  000  038  112  148  127  078  036  013 


2  160  150  019  023  095  076  010  015  087  134  117  070 


3  223  041  029  097  027  010  075  066  006  022  095  128 


4  221  003  102  021  027  076  015  018  076  044  000  046 


5  163  082  054  019  074  003  042  056  001  042  070  015 


6  091  164  000  080  007  046  041  004  059  023  009  066 


7  039  165  063  036  030  047  005  058  006  033  048  000 


8  012  106  152  003  066  000  058  004  039  029  008  055 


9  002  045  157  084  011  047  014  035  026  014  046  000 


10 

000 

012 

094 

161 

023 

040 

015 

040 

010 

043 

001 

047 

11 

000 

001 

034 

140 

124 

000 

049 

000 

050 

000 

044 

003 

12 

000 

000 

006 

068 

164 

071 

008 

037 

007 

043 

005 

033 

000  000  017  107  159  028  022  017  023  027  017 


14  000  000  000  001  036  141  133  006  030  004  038  012 


15 

000 

000 

000 

000 

004 

061 

163 

100 

000 

028 

000 

046 

16 

000 

000 

000 

001 

000 

010 

088 

169 

068 

002 

020 

005 

17 

000 

000 

000 

001 

002 

000 

018 

113 

164 

043 

006 

Oil 

18 

000 

000 

000 

000 

001 

003 

000 

028 

133 

151 

026 

008 

19 

000 

000 

000 

000 

000 

003 

004 

000 

040 

147 

135 

014 

24  000  000  000  000  000  000  000  000  001  008  021  013 


Values  in  the  table  are  (FCF  x  1000) 


TABLE  4  (cont.) 

Franck-Condon  Factors  for  PbO  Small  B  to  X  Transitions 


V" 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

V’ 

0 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

1 

004 

001 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

2 

031 

010 

003 

001 

000 

000 

000 

000 

000 

000 

000 

000 

3 

100 

053 

020 

006 

000 

000 

000 

000 

000 

000 

000 

000 

4 

113 

118 

076 

033 

002 

000 

000 

000 

000 

000 

000 

000 

5 

012 

084 

123 

095 

016 

004 

001 

000 

000 

000 

000 

000 

6 

040 

000 

053 

117 

061 

022 

005 

001 

000 

000 

000 

000 

7 

044 

059 

005 

028 

120 

074 

028 

006 

001 

000 

000 

000 

8 

011 

021 

065 

018 

091 

126 

085 

032 

007 

001 

000 

000 

9 

045 

027 

006 

060 

003 

077 

130 

093 

036 

008 

001 

000 

10 

007 

028 

041 

000 

044 

000 

066 

132 

100 

039 

008 

001 

11 

036 

020 

012 

048 

038 

051 

000 

057 

135 

107 

041 

007 

12 

015 

021 

032 

003 

008 

028 

056 

002 

052 

138 

112 

041 

13 

018 

027 

009 

039 

043 

014 

019 

058 

004 

048 

142 

115 

14 

027 

007 

035 

002 

001 

037 

021 

013 

058 

005 

048 

147 

15 

003 

032 

001 

038 

038 

005 

031 

025 

010 

058 

006 

051 

16 

046 

000 

032 

000 

001 

034 

009 

026 

028 

008 

057 

005 

17 

015 

041 

001 

028 

034 

004 

029 

012 

022 

029 

007 

055 

18 

004 

027 

034 

004 

007 

030 

007 

025 

015 

019 

030 

006 

19 

009 

000 

038 

026 

017 

012 

026 

009 

022 

016 

018 

029 

20 

007 

008 

001 

046 

Oil 

012 

016 

023 

Oil 

019 

017 

018 

21 

100 

004 

006 

004 

015 

012 

008 

019 

020 

012 

017 

017 

22 

160 

084 

002 

004 

055 

012 

013 

005 

022 

019 

012 

016 

23 

079 

155 

070 

001 

015 

058 

009 

012 

003 

025 

019 

Oil 

24 

004 

084 

147 

057 

001 

021 

060 

008 

Oil 

002 

026 

019 

Values  in  the  table  are  (FCF  x  1000) 
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TABLE  5 


Franck-Condon  Factors  for  PbO  Large  A  to  X  Transitions 


V" 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

.  *  . 

V' 

0 

026 

097 

179 

218 

197 

141 

082 

039 

015 

005 

001 

000 

• 

1 

091 

180 

127 

022 

009 

083 

148 

149 

104 

054 

022 

007 

w-  -■  - 

2 

165 

129 

004 

051 

108 

049 

000 

048 

122 

139 

101 

053 

• 

3 

204 

028 

046 

095 

Oil 

029 

090 

047 

000 

048 

119 

129 

•;  v- 

4 

190 

003 

102 

013 

040 

074 

005 

035 

082 

030 

003 

065 

‘  . * 

5 

143 

058 

058 

020 

074 

001 

051 

052 

000 

053 

070 

Oil 

6 

090 

116 

003 

077 

Oil 

043 

044 

004 

062 

025 

010 

069 

c 

7 

049 

131 

017 

061 

015 

058 

001 

058 

012 

026 

057 

003 

8 

024 

108 

070 

010 

064 

006 

049 

019 

024 

046 

000 

052 

9 

010 

072 

106 

006 

056 

016 

044 

009 

049 

000 

050 

013 

10 

004 

041 

105 

047 

012 

058 

002 

051 

004 

042 

014 

024 

*  .*.• 

11 

002 

020 

080 

087 

003 

050 

019 

030 

021 

029 

013 

039 

12 

001 

009 

051 

098 

036 

012 

054 

000 

048 

000 

044 

000 

*• 

• 

13 

000 

004 

028 

081 

075 

002 

043 

024 

017 

033 

010 

032 

14 

000 

001 

013 

055 

090 

032 

009 

052 

001 

039 

008 

031 

%... 

15 

000 

000 

006 

032 

079 

068 

003 

035 

031 

006 

039 

000 

16 

000 

000 

002 

016 

056 

084 

031 

005 

048 

007 

026 

021 

17 

000 

000 

001 

007 

033 

075 

064 

004 

026 

037 

001 

038 

18 

000 

000 

000 

003 

017 

054 

078 

033 

002 

043 

015 

014 

yo« 

A 

19 

000 

000 

000 

001 

008 

033 

070 

062 

007 

018 

040 

001 

ii> 

w 

20 

000 

000 

000 

000 

003 

018 

051 

073 

035 

000 

036 

023 

L 

21 

000 

000 

000 

000 

001 

008 

032 

065 

061 

Oil 

Oil 

040 

22 

000 

000 

000 

000 

000 

004 

017 

047 

069 

039 

000 

028 

23 

000 

000 

000 

000 

000 

001 

008 

029 

060 

060 

016 

005 

•VO 

24 

000 

000 

000 

000 

000 

001 

004 

016 

043 

065 

042 

002 

Values  in  the  table  are  (FCF  x  1000)  • 


4 


TABLE  5 


( cont . ) 


Franck-Condon  Factors  for  PbO  Large  A  to  X  Transitions 


V" 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

V' 

0 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

1 

002 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

2 

021 

007 

002 

000 

000 

000 

000 

000 

000 

000 

000 

000 

3 

088 

044 

017 

005 

000 

000 

000 

000 

000 

000 

000 

000 

4 

123 

116 

072 

032 

003 

001 

000 

000 

000 

000 

000 

000 

5 

016 

089 

125 

099 

020 

006 

001 

000 

000 

000 

000 

000 

6 

046 

000 

044 

112 

076 

034 

Oil 

003 

000 

000 

000 

000 

7 

036 

069 

016 

010 

124 

099 

051 

018 

005 

001 

000 

000 

8 

029 

005 

063 

046 

046 

115 

117 

070 

028 

008 

002 

000 

9 

022 

005 

002 

035 

012 

016 

093 

125 

089 

040 

012 

003 

10 

039 

001 

051 

022 

065 

035 

001 

065 

125 

107 

054 

018 

11 

002 

047 

007 

028 

000 

047 

056 

002 

037 

115 

120 

069 

12 

043 

005 

031 

028 

052 

010 

023 

065 

015 

016 

099 

130 

13 

013 

024 

025 

009 

000 

041 

028 

006 

062 

033 

003 

079 

14 

009 

033 

004 

039 

041 

Oil 

022 

044 

000 

048 

049 

000 

15 

039 

000 

038 

001 

008 

025 

028 

006 

049 

006 

030 

060 

16 

013 

026 

010 

025 

020 

024 

009 

040 

000 

044 

019 

014 

17 

003 

030 

007 

027 

029 

005 

036 

000 

041 

005 

031 

033 

18 

031 

002 

034 

000 

000 

035 

000 

036 

003 

032 

016 

016 

19 

030 

013 

015 

022 

027 

006 

028 

007 

026 

014 

018 

029 

20 

004 

034 

001 

028 

021 

012 

018 

014 

019 

012 

026 

007 

21 

005 

019 

024 

003 

000 

030 

002 

028 

003 

029 

003 

033 

22 

030 

000 

030 

008 

021 

005 

027 

001 

030 

000 

031 

000 

23 

037 

012 

009 

029 

026 

008 

016 

017 

009 

023 

005 

026 

24 

020 

034 

001 

022 

004 

028 

000 

025 

006 

019 

012 

015 

Values  in  the  table  are  (FCF  x  1000) 


C 


C 
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c 

[ 


TABLE  6 


Franck-Condon  Factors 


for  PbO  Large  B 


to  X  Transitions 


e 


c 


V" 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

V' 

0 

049 

160 

244 

235 

165 

090 

039 

014 

004 

001 

000 

000 

1 

126 

196 

073 

001 

079 

166 

168 

112 

053 

019 

005 

001 

2 

176 

090 

006 

112 

086 

003 

045 

139 

159 

109 

051 

018 

3 

178 

009 

085 

069 

003 

090 

074 

002 

047 

136 

148 

095 

■m 

150 

008 

101 

000 

080 

044 

009 

089 

054 

000 

065 

140 

5 

111 

045 

049 

037 

063 

005 

080 

018 

028 

089 

029 

008 

6 

076 

075 

007 

079 

005 

064 

024 

028 

068 

001 

056 

075 

7 

049 

084 

002 

067 

014 

059 

007 

066 

001 

059 

037 

008 

8 

030 

076 

019 

031 

052 

012 

054 

014 

040 

036 

012 

067 

9 

018 

061 

035 

006 

064 

002 

058 

006 

054 

002 

060 

004 

10 

Oil 

046 

043 

000 

049 

026 

022 

045 

Oil 

043 

016 

031 

11 

007 

033 

043 

005 

026 

046 

001 

057 

004 

047 

005 

049 

12 

004 

023 

039 

012 

009 

048 

006 

035 

033 

013 

041 

008 

13 

003 

017 

032 

016 

002 

039 

020 

010 

051 

001 

046 

007 

14 

002 

012 

026 

018 

000 

026 

030 

000 

044 

018 

019 

037 

15 

001 

009 

021 

018 

001 

015 

033 

002 

026 

037 

001 

047 

16 

001 

006 

017 

017 

002 

008 

029 

009 

010 

043 

005 

029 

17 

001 

005 

013 

015 

004 

003 

023 

014 

002 

035 

018 

008 

18 

001 

004 

Oil 

013 

005 

001 

017 

017 

000 

024 

029 

000 

19 

000 

003 

009 

012 

005 

000 

012 

017 

001 

013 

031 

003 

20 

000 

002 

007 

010 

005 

000 

009 

015 

003 

006 

028 

010 

21 

000 

002 

006 

009 

005 

000 

006 

013 

005 

002 

022 

016 

22 

000 

002 

005 

007 

005 

000 

004 

Oil 

006 

000 

015 

018 

23 

000 

001 

004 

006 

004 

000 

003 

009 

006 

000 

010 

018 

24 

000 

001 

004 

005 

004 

000 

002 

007 

006 

000 

007 

016 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  the  Simpson’s  Rule  Approximation 


C 


c 


44 


< 


V' 

0 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

1  -•! 

m 

1 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

2 

005 

001 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

3 

043 

015 

004 

001 

000 

000 

000 

000 

000 

000 

000 

000 

4 

133 

079 

033 

Oil 

000 

000 

000 

000 

000 

000 

000 

000 

J 

5 

090 

142 

115 

061 

007 

001 

000 

000 

000 

000 

000 

000 

6 

007 

030 

115 

137 

042 

014 

003 

001 

000 

000 

000 

000 

7 

077 

046 

001 

065 

123 

068 

026 

007 

001 

000 

000 

000 

V 

8 

006 

036 

077 

015 

105 

139 

098 

044 

014 

003 

000 

000 

c 

9 

046 

043 

002 

067 

000 

059 

135 

125 

068 

025 

006 

001 

'  m  *  «  * 

10 

039 

007 

062 

010 

077 

016 

018 

107 

142 

096 

041 

012 

11 

001 

056 

006 

039 

001 

066 

051 

000 

064 

141 

122 

063 

12 

038 

017 

028 

038 

059 

Oil 

030 

075 

012 

023 

118 

143 

13 

039 

005 

046 

001 

005 

038 

041 

002 

069 

043 

002 

081 

9 

14 

006 

039 

006 

038 

029 

033 

007 

058 

007 

038 

070 

006 

c 

15 

005 

035 

009 

034 

043 

002 

050 

002 

042 

034 

008 

073 

16 

029 

008 

037 

003 

004 

041 

009 

034 

025 

013 

056 

001 

17 

044 

001 

036 

009 

012 

029 

012 

036 

006 

045 

000 

051 

18 

038 

018 

013 

034 

036 

001 

040 

000 

044 

002 

03° 

015 

19 

022 

035 

000 

038 

024 

014 

019 

020 

020 

022 

024 

014 

--v 

*  .  *  -V 

!• 

20 

008 

039 

005 

022 

003 

034 

000 

035 

001 

039 

001 

041 

> 

21 

001 

032 

018 

006 

003 

026 

014 

014 

022 

010 

030 

008 

22 

000 

021 

028 

000 

019 

007 

032 

000 

030 

002 

031 

007 

23 

002 

Oil 

030 

004 

033 

000 

029 

012 

013 

022 

005 

031 

24 

005 

005 

027 

Oil 

035 

008 

014 

029 

000 

028 

003 

024 

0  006 


1  024 


2  054 


3  089 


4  117 


5  132 


6  130 


037  105  184  223  199  135  071  029  009  002  000 


098  154  105  015  016  106  173  157  094  040  013 


134  089  002  047  104  045  001  072  154  150  091 


118  013  035  084  014  024  091  045  001  074  149 


068  005  076  020  022  072  011  030  087  029  008 


021  042  049  004  063  011  030  061  001  049  077 


000  068  007  043  029  014  054  001  049  040  003 


:  7 

117 

008 

059 

004 

053 

e*ei 

049 

008 

033 

034 

006 

061  | 

8 

097 

032 

030 

029 

025 

021 

031 

010 

042 

001 

049 

008  1 

9  075  056  006  048  001  043  002  039  006  034  017  020 


10  056  073  000  044  006  033  009  031  008  034  004  041 


11  039  078  013  024  027  009  031  005  032  004  032  007 


12  026  075  034  006  039  000  033  004  028  007  027  007 


13  017  064  054  000  033  012  015  022  006  028  002  030 


14  010  050  067  011  015  028  001  030  002  025  007  021 


15  006  036  069  031  002  032  005  018  017  006  024  001 


16  003  024  061  051  002  020  019  003  026  001  022  007 


17  001  014  049  063  017  006  027  002  018  013  006  022 


18  001  008  035  065  038  000  021  012  004  023  001  019 


21 

000 

001 

007 

030 

062 

049 

006 

008 

017 

000 

017 

010 

22 

000 

000 

003 

018 

051 

062 

023 

000 

016 

006 

006 

019 

23 

000 

000 

001 

009 

036 

063 

045 

004 

007 

013 

000 

016 

v-vl 

24 

000 

000 

000 

004 

022 

054 

061 

021 

000 

012 

004 

006 

TABLE  7  (cont.) 


Franck-Condon  Factors  for  PbO  Large  C  to  X  Transitions 


V” 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

V’ 

0 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

1 

003 

001 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

2 

039 

013 

003 

001 

000 

000 

000 

000 

000 

000 

000 

000 

3 

137 

079 

032 

010 

000 

000 

000 

000 

000 

000 

000 

000 

4 

092 

148 

121 

063 

007 

001 

000 

000 

000 

000 

000 

000 

5 

010 

025 

114 

143 

046 

016 

004 

001 

000 

000 

000 

000 

6 

070 

054 

000 

055 

130 

077 

032 

009 

002 

000 

000 

000 

7 

013 

024 

078 

024 

091 

140 

109 

055 

020 

005 

001 

000 

8 

031 

050 

000 

055 

003 

039 

121 

131 

083 

036 

012 

003 

9 

040 

001 

054 

021 

074 

032 

004 

078 

134 

111 

059 

022 

10 

000 

044 

012 

023 

000 

048 

065 

005 

031 

112 

129 

086 

11 

028 

018 

016 

040 

051 

022 

012 

071 

033 

003 

071 

128 

12 

030 

003 

038 

000 

012 

023 

049 

000 

047 

063 

005 

029 

13 

002 

031 

005 

029 

018 

037 

001 

051 

018 

014 

069 

030 

14 

010 

020 

010 

025 

035 

000 

043 

007 

028 

044 

000 

049 

15 

028 

000 

030 

000 

001 

033 

009 

024 

030 

004 

052 

013 

16 

016 

013 

012 

017 

017 

018 

012 

029 

003 

043 

003 

034 

17 

000 

025 

000 

025 

027 

001 

031 

000 

036 

003 

032 

022 

18 

007 

012 

015 

006 

004 

023 

005 

025 

009 

021 

020 

010 

19 

020 

000 

022 

002 

005 

019 

008 

020 

007 

025 

003 

034 

20 

016 

007 

009 

016 

021 

000 

024 

000 

027 

000 

030 

002 

21 

004 

018 

000 

018 

014 

009 

Oil 

014 

009 

017 

010 

016 

22 

001 

014 

008 

006 

001 

020 

000 

021 

002 

021 

003 

024 

23 

009 

003 

016 

000 

005 

010 

012 

005 

018 

002 

021 

001 

24 

017 

001 

Oil 

008 

016 

000 

017 

002 

015 

006 

013 

010 

Values  in  the  table  are  (FCF  x  1000) 


TABLE  8 


Franck-Condon  Factors  for  PbO  C  Prime  to  X  Transitions 


V" 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

V' 

0 

002 

014 

051 

110 

171 

200 

183 

134 

078 

037 

014 

004 

1 

010 

052 

117 

138 

082 

Oil 

012 

085 

149 

151 

105 

055 

2 

027 

097 

122 

050 

000 

052 

095 

047 

000 

043 

120 

144 

3 

052 

119 

065 

000 

053 

072 

010 

021 

082 

054 

001 

036 

4 

080 

105 

010 

033 

068 

007 

029 

066 

012 

019 

078 

047 

5 

103 

065 

003 

067 

019 

018 

059 

006 

029 

059 

006 

028 

6 

116 

025 

032 

050 

001 

054 

013 

021 

050 

001 

040 

047 

7 

117 

003 

058 

014 

030 

035 

005 

049 

003 

034 

035 

002 

8 

109 

002 

059 

000 

050 

003 

037 

017 

017 

038 

001 

046 

9 

095 

017 

040 

014 

036 

007 

039 

001 

042 

002 

035 

018 

10 

079 

037 

017 

036 

010 

031 

012 

024 

019 

015 

030 

004 

11 

062 

055 

002 

045 

000 

039 

000 

037 

000 

036 

001 

034 

12 

047 

066 

001 

038 

Oil 

024 

015 

019 

016 

019 

014 

023 

13 

035 

069 

013 

020 

028 

005 

032 

001 

032 

000 

032 

000 

14 

024 

066 

029 

005 

037 

001 

030 

006 

022 

Oil 

018 

013 

15 

016 

057 

045 

000 

032 

012 

013 

022 

004 

027 

001 

028 

16 

010 

046 

055 

006 

018 

026 

001 

029 

002 

023 

008 

017 

17 

006 

035 

058 

019 

005 

032 

003 

019 

015 

007 

023 

001 

18 

004 

025 

054 

034 

000 

026 

014 

005 

025 

000 

022 

005 

19 

002 

017 

047 

045 

005 

014 

025 

000 

022 

009 

009 

019 

20 

001 

Oil 

038 

051 

016 

003 

027 

006 

010 

020 

000 

021 

21 

001 

007 

029 

051 

029 

000 

020 

017 

001 

022 

005 

010 

22 

000 

004 

021 

046 

040 

005 

009 

024 

002 

013 

016 

001 

23 

000 

002 

014 

038 

047 

016 

002 

022 

010 

003 

020 

003 

24 

000 

001 

009 

030 

047 

028 

001 

014 

019 

000 

015 

012 

Values  in  the  table  are  (FCF  x  1000) 


TABLE  8  (cont.) 

Franck-Condon  Factors  for  PbO  C  Prime  to  I  Transitions 


V"  12  13  14  15  16  17  18  19  20  21  22  23 


2 

109 

059 

025 

008 

000 

000 

000 

000 

000 

000 

000 

000 

3 

112 

136 

102 

054 

007 

002 

000 

000 

000 

000 

000 

000 

4 

000 

044 

116 

130 

045 

016 

005 

001 

000 

000 

000 

000 

5 

076 

031 

002 

062 

121 

075 

033 

Oil 

003 

001 

000 

000 

6 

000 

044 

069 

014 

087 

131 

106 

056 

021 

006 

001 

000 

c 

7 

052 

028 

005 

062 

001 

038 

113 

127 

084 

038 

012 

003 

8 

015 

016 

055 

008 

070 

026 

005 

073 

128 

111 

060 

023 

9 

013 

044 

001 

038 

000 

051 

058 

004 

030 

108 

127 

086 

10 

041 

001 

035 

025 

052 

014 

016 

068 

029 

00? 

071 

128 

•VV 

11 

007 

024 

025 

006 

003 

032 

041 

000 

048 

058 

004 

032 

12 

008 

031 

001 

038 

030 

026 

005 

051 

013 

018 

067 

025 

A 

13 

032 

002 

029 

009 

025 

006 

042 

002 

033 

038 

000 

052 

Mm 

14 

018 

012 

022 

007 

001 

036 

002 

032 

022 

009 

050 

007 

15 

000 

029 

000 

030 

029 

006 

023 

020 

009 

039 

000 

039 

.*  \\ 

16 

012 

013 

015 

014 

017 

010 

025 

004 

034 

000 

037 

013 

-■  ~ 

17 

025 

000 

026 

000 

000 

029 

000 

031 

002 

029 

Oil 

018 

18 

015 

012 

010 

016 

019 

007 

019 

010 

017 

017 

Oil 

029 

*  19 

001 

023 

000 

023 

022 

003 

022 

003 

025 

002 

029 

000 

-*  V  * 

*  20 

004 

014 

Oil 

008 

004 

021 

002 

023 

001 

025 

003 

025 

C 

21 

016 

002 

021 

000 

003 

017 

007 

014 

010 

014 

on 

016 

22 

020 

003 

013 

Oil 

018 

001 

021 

000 

022 

000 

024 

001 

V*/' 

23 

011 

013 

002 

019 

017 

005 

012 

Oil 

007 

016 

005 

019 

24 

002 

018 

002 

013 

004 

017 

000 

019 

001 

018 

004 

017 

2 

196 

032 

064 

128 

007 

071 

198 

178 

089 

029 

007 

001 

m  .  i 

*  .* 

3 

174 

002 

112 

008 

079 

094 

000 

096 

192 

149 

068 

021 

»  .*  .  J 

4 

133 

038 

060 

027 

081 

001 

102 

052 

010 

130 

178 

117 

•y*  *-#i 

5 

092 

079 

008 

074 

008 

068 

034 

030 

104 

013 

042 

151 

6 

059 

097 

003 

060 

013 

059 

008 

079 

001 

076 

071 

000 

v  V 

c 

7 

036 

093 

028 

021 

049 

008 

056 

012 

052 

042 

016 

096 

8 

021 

076 

056 

001 

053 

005 

047 

Oil 

053 

006 

074 

004 

9 

012 

057 

072 

008 

028 

032 

009 

046 

004 

053 

010 

046 

10 

006 

039 

074 

030 

005 

042 

002 

039 

012 

036 

015 

048 

o-::: 

11 

003 

026 

066 

053 

001 

029 

020 

010 

038 

002 

047 

001 

>> 

12 

002 

016 

053 

066 

016 

008 

032 

000 

033 

Oil 

025 

021 

• 

13 

001 

009 

038 

068 

040 

000 

025 

013 

010 

032 

000 

039 

14 

000 

004 

024 

060 

060 

Oil 

008 

024 

000 

028 

Oil 

017 

15 

000 

002 

014 

046 

068 

034 

000 

019 

008 

009 

029 

000 

16 

000 

001 

007 

030 

064 

056 

009 

006 

017 

000 

025 

012 

17 

000 

000 

003 

017 

050 

068 

032 

000 

013 

005 

009 

027 

88 

18 

000 

000 

001 

007 

032 

065 

056 

Oil 

003 

010 

000 

024 

• 

19 

000 

000 

000 

002 

017 

050 

068 

035 

001 

007 

002 

010 

20 

000 

000 

000 

000 

006 

031 

064 

059 

016 

000 

004 

001 

21 

000 

000 

poo 

000 

001 

013 

046 

068 

042 

005 

002 

000 

22 

000 

000 

000 

001 

000 

003 

023 

057 

061 

024 

000 

001 

23 

000 

000 

000 

001 

001 

000 

006 

033 

060 

046 

010 

000 

24 

000 

000 

000 

001 

002 

002 

000 

Oil 

040 

054 

028 

003 

h  V’ 

Values  in  the  table  are  (FCF 


TABLE  9  (cont.) 

Franck-Condon  Factors  for  PbO  Large  D  to  X  Transitions 
V"  12  13  14  15  16  17  18  19  20  21  22  23 


I 

0  000  000  000  000  000  000  000  000  000  000  000  000 


1  000  000  000  000  000  000  000  000  000  000  000  000 


2  000  000  000  000  000  000  000  000  000  000  000  000 


3  005  001  000  000  000  000  000  000  000  000  000  000 


4  050  016  004  001  000  000  000  000  000  000  000  000 


5  155  091  037  012  001  000  000  000  000  000  000  000 


6  078  154  131  071  009  002  001  000  000  000  000  000 


7  029  013  104  146  057  023  007  002  000  000  000  000 


8  056  079  005  035  134  093  047  019  006  002  000  000 


9 

046 

006 

080 

049 

056 

122 

122 

080 

039 

016 

005 

001 

10 

004 

068 

Oil 

032 

023 

007 

072 

121 

110 

070 

034 

014 

11 

052 

011 

036 

052 

056 

065 

008 

017 

082 

117 

100 

062 

12  028  017  046  001  023  010  067  046  001  028  087  112 


13  000  046  000  047  021  055  005  026  066  030  000  036 


14  022  015  026  024  046  000  042  038  000  039  060  019 


15 _ 032  001  037  001  001  038  023  007  049  020  004  047 


19  027  001  017  002  008  022  003  032  000  027  016  005 


20  025  015  003  013  028  001  022  007  020  017  003  031 


21  013  030  002  010  023  010  014  003  027  001  022  008 


22  003  029  019  001  005  028  000  015  005  021  009  004 


23  000  017  035  006  000  023  013  007  002  025  002  015 


24  000  006  033  024  003  007  030  000  008  007  021  004 


Values  in  the  table  are  (FCF 


1 

003 

017 

052 

095 

118 

098 

042 

002 

020 

080 

128 

132 

'  V ' 

2 

Oil 

051 

098 

104 

055 

005 

014 

065 

075 

029 

000 

033 

A 

ESS 

w 

3 

033 

094 

104 

043 

000 

028 

062 

033 

000 

032 

068 

043 

5 

4 

074 

117 

054 

000 

031 

054 

016 

005 

045 

040 

002 

020 

5 

128 

090 

005 

026 

051 

012 

008 

042 

021 

002 

037 

038 

r 

6 

179 

031 

Oil 

054 

015 

007 

038 

014 

005 

036 

018 

002 

a 

7 

201 

000 

043 

029 

003 

034 

014 

005 

031 

010 

006 

033 

► .  *• 

8 

178 

045 

035 

001 

034 

017 

002 

029 

008 

007 

029 

005 

• 

9 

117 

137 

001 

009 

035 

001 

020 

012 

006 

022 

004 

012 

9  , 

10 

049 

181 

035 

010 

008 

030 

008 

001 

025 

002 

008 

021 

.*  'i 

11 

008 

125 

131 

001 

000 

042 

004 

008 

Oil 

012 

010 

001 

11:; 

HR 

• 

12 

000 

036 

155 

048 

000 

020 

040 

000 

000 

029 

000 

004 

iu. 

13 

005 

000 

075 

108 

009 

005 

053 

015 

000 

016 

023 

001 

14 

005 

Oil 

005 

081 

039 

002 

026 

048 

003 

003 

038 

008 

*.**“.* 

e 

15 

001 

017 

012 

on 

041 

004 

004 

038 

019 

001 

018 

030 

1 

TABLE  10  (cont.) 

Franck-Condon  Factors  for  PbO  Large  E  to  X  Transitions 
V"  12  13  14  15  16  17  18  19  20  21  22  23 


V11 

0 

017 

007 

002 

001 

000 

000 

000 

000 

000 

000 

000 

000 

1 

102 

062 

031 

013 

001 

000 

000 

000 

000 

000 

000 

000 

2 

092 

123 

108 

072 

017 

007 

002 

001 

000 

000 

000 

000 

3 

003 

017 

072 

110 

075 

043 

021 

009 

003 

001 

000 

000 

4 

060 

045 

005 

Oil 

098 

100 

076 

046 

024 

Oil 

004 

002 

5 

003 

017 

056 

043 

008 

051 

089 

095 

075 

049 

027 

013 

6 

035 

034 

002 

018 

040 

006 

006 

044 

080 

089 

075 

052 

7 

013 

003 

035 

029 

019 

049 

037 

006 

004 

036 

070 

084 

8 

009 

031 

008 

006 

024 

000 

018 

045 

035 

007 

002 

027 

9 

023 

002 

015 

027 

010 

033 

020 

000 

016 

041 

035 

010 

10 

000 

014 

019 

000 

024 

001 

010 

032 

018 

000 

013 

035 

11 

018 

009 

000 

021 

000 

021 

017 

001 

Oil 

026 

019 

001 

12 

018 

003 

009 

Oil 

013 

Oil 

003 

015 

020 

001 

005 

026 

13 

002 

026 

000 

000 

001 

001 

025 

003 

000 

021 

Oil 

006 

14 

000 

023 

014 

001 

015 

003 

010 

Oil 

005 

016 

003 

003 

15 

001 

006 

033 

002 

030 

002 

000 

028 

000 

001 

025 

001 

Values  in  the  table  are  (FCF  x  1000) 


53 


4' 


TABLE  11 

Franck-Condon  Factors  for  PbO  A  to  Small  A  Transitions 


2  3 


7  8 


10  11 


V' 

0 

918 

070 

010 

002 

000 

000 

000 

000 

000 

000 

000 

000 

• 

1 

082 

765 

120 

026 

006 

001 

000 

000 

000 

000 

000 

000 

2 

000 

164 

623 

154 

043 

012 

003 

001 

000 

000 

000 

000 

•  - 

1 

3 

000 

000 

244 

493 

174 

060 

020 

006 

002 

000 

000 

000 

V' 

4 

000 

000 

002 

320 

378 

182 

075 

029 

010 

003 

001 

000 

5 

000 

000 

000 

004 

390 

277 

180 

086 

038 

015 

005 

002 

6 

000 

000 

000 

000 

009 

451 

193 

171 

094 

048 

021 

008 

C' 

7 

000 

000 

000 

000 

001 

016 

502 

125 

157 

098 

056 

027 

8 

000 

000 

000 

000 

000 

001 

026 

541 

071 

140 

099 

063 

9 

000 

000 

000 

000 

000 

000 

001 

041 

567 

033 

122 

097 

10 

000 

000 

000 

000 

000 

000 

000 

000 

061 

576 

010 

103 

1 

11 

000 

000 

000 

000 

000 

000 

000 

000 

000 

087 

567 

000 

12 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

118 

542 

*  **, 

A 

13 

000 

000 

000 

000 

000 

000 

000 

000 

000 

001 

000 

155 

w 

14 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

001 

002 

P 

15 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

001 

*  *  •  * 

16 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

17 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

,*/,v 

18 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

19 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

t 

% 

20 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

c 

21 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

•yv 

22 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

23 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

24 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

Values  in  the  table  are  (FCF 


V’ 

0 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

• 

1 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

-  —  < 

2 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

•  ,  J 

3 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

A 

4 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

5 

001 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

6 

003 

001 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

c 

7 

012 

005 

002 

001 

000 

000 

000 

000 

000 

000 

000 

000 

8 

033 

016 

007 

002 

000 

000 

000 

000 

000 

000 

000 

000 

iiiiiiiigiB 

9 

068 

039 

020 

009 

001 

000 

000 

000 

000 

000 

000 

000 

10 

093 

071 

045 

024 

004 

002 

001 

000 

000 

000 

000 

000 

11 

085 

088 

073 

050 

013 

006 

002 

001 

000 

000 

000 

000 

* .  * 

12 

004 

067 

082 

073 

032 

016 

007 

003 

001 

000 

000 

000 

• 

13 

501 

018 

051 

076 

057 

036 

019 

009 

004 

001 

000 

000 

£ 

14 

194 

447 

042 

035 

070 

059 

040 

023 

Oil 

005 

002 

001 

15 

004 

235 

382 

073 

061 

067 

060 

044 

026 

014 

006 

002 

-  **.y 

16 

002 

009 

274 

311 

010 

054 

063 

061 

047 

030 

016 

008 

17 

000 

003 

017 

309 

148 

003 

045 

060 

060 

050 

033 

019 

*•* 

18 

000 

001 

005 

028 

166 

184 

000 

037 

056 

059 

052 

037 

♦  •  ***  i 

I  A 

19 

000 

000 

001 

008 

353 

101 

214 

003 

029 

051 

057 

053 

20 

000 

000 

000 

002 

061 

356 

048 

234 

014 

020 

047 

055 

c 

21 

000 

000 

000 

000 

018 

081 

343 

014 

240 

031 

012 

043 

22 

000 

000 

000 

000 

005 

026 

104 

314 

000 

230 

055 

005 

: 

23 

000 

000 

000 

000 

000 

008 

037 

125 

270 

008 

204 

084 

24 

000 

000 

000 

000 

000 

001 

012 

051 

143 

215 

035 

164 

•  v%*« 

V' 


0 

526 

330 

114 

026 

004 

000 

000 

000 

000 

000 

000 

000 

1 

358 

081 

294 

190 

063 

012 

001 

000 

000 

000 

000 

000 

2 

102 

358 

000 

194 

218 

099 

025 

003 

000 

000 

000 

000 

3 

013 

195 

276 

017 

112 

213 

127 

040 

007 

000 

000 

000 

4 

001 

034 

258 

201 

045 

058 

191 

145 

055 

Oil 

001 

000 

5 

000 

002 

056 

297 

150 

064 

027 

163 

153 

069 

017 

002 

6 

000 

000 

002 

072 

325 

122 

071 

012 

136 

154 

080 

023 

7 

000 

000 

000 

002 

082 

348 

113 

068 

004 

113 

150 

087 

8 

000 

000 

000 

000 

001 

082 

369 

118 

060 

001 

096 

142 

9 

000 

000 

000 

000 

001 

000 

073 

388 

137 

047 

000 

084 

10 

000 

000 

000 

000 

000 

001 

001 

055 

399 

170 

031 

000 

11 

000 

000 

000 

000 

000 

000 

002 

005 

032 

396 

217 

016 

12 

000 

000 

000 

000 

000 

000 

000 

002 

014 

Oil 

369 

275 

13 

000 

000 

000 

000 

000 

000 

000 

000 

001 

025 

000 

313 

14 

000 

000 

000 

000 

000 

000 

000 

000 

001 

000 

036 

010 

15 

000 

000 

000 

000 

000 

000 

000 

000 

000 

002 

000 

039 

16 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

002 

005 

17 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

002 

18 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

002 

19 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

20 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

21 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

22 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

23 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

24 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

Values  in  the  table  are  (FCF  x  1000) 


C 


c 
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TABLE  12  (cont.) 

Franck-Condon  Factors  for  PbO  D  to  Small  A  Transitions 
V"  12  13  14  15  16  17  18  19  20  21  22  23 


0  000 


1  000 


2  000 


3  000 


4  000 


5  000 


6  003 


7  028 


8  092 


9  134 


10  078 


11  000 


12  004 


13  337 


14  227 


15  046 


16  031 


17  017 


18  000 


19  003 


20  001 


21  000 


22  000 


23  000 


24  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


005  000 


032  006 


093  035 


126  093 


076  117 


001  078 


000  004 


387  009 


125  406 


104  037 


012  159 


032  000 


001  039 


003  011 


002  001 


000  004 


000  001 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


007  000 


037  001 


091  009 


106  038 


085  088 


010  076 


030  121 


372  044 


000  068 


174  156 


018  128 


027  043 


029  125 


001  006 


004  028 


005  036 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


000  000 


002  000 


014  003 


035  017 


101  029 


013  118 


243  001 


072  264 


000  035 


004  027 


079  005 


094  011 


001  082 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


000  000  000 


001  000  000 


004  001  000 


020  005  002 


021  025  005 


147  011  035 


002  188  002 


234  012  232 


002  144  015 


082  012  042 


001  091  047 


000  042  036 


Values  in  the  table  are  (FCF 


V 


c 


( 


c 


Molecular  Bromine 

Three  isotopic  forms  of  molecular  Bromine  were  studied, 
the  (79,79)  isotope,  the  (81,81)  isotope,  and  the  (79,81) 
mixture  of  both  isotopic  states.  Five  electronic  states  were 
calculated  for  the  first  isotope  (A',  B,  D',  E,  X)  and  two 
states  each  (B,  X)  for  the  (81,81)  isotope  and  the  mixed 
isotope.  The  data  used  to  calculate  the  bromine  potentials 
is  given  in  Table  13.  The  number  next  to  the  state  name  is 
the  bibliographic  entry  for  the  journal  article  from  which 
the  molecular  constants  were  drawn.  The  dissociation  energy 
for  the  B  and  X  states  of  bromine  are  difficult  to  calculate 
and  not  widely  available  in  the  open  literature.  Because  of 
this,  three  spectroscopic  dissociation  energies  (79:B,  81:B, 
79,81:X)  had  to  be  calculated  from  the  chemical  dissociation 
energies  which  were  published  in  the  article  by  Barrow  and 
others  (2).  The  formula  used  to  make  the  conversion  was: 

De  =  Do  +  1/2ue  ~  1/4u)eXe  +  1/8weYe  (22) 

where  De  is  the  desired  spectroscopic  dissociation  energy,  D0 
is  the  chemical  dissociation  energy  and  U)e,  weXe  ,  we  Ye  are 
the  molecular  constants  described  in  Section  II.  Three 
dissociation  energies  (79:X,  81:X,  79,81:B)  had  to  be 
extrapolated  from  available  data  because  no  specific  energy 
data  was  available.  The  method  of  extrapolation  is  described 
in  Appendix  C.  The  remaining  three  dissociation  energies 
were  available  in  the  literature. 
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TABLE  13 


Data  for  Molecular  Bromine 
Part  I.  (79,79) 


STATE 

A’  (26) 

B  (2) 

D'  (26) 

E  (3) 

X  (2) 

T 

Ae 

13220 

15902.47 

48930 

49779.06 

0 

165.172 

167.6066 

150.863 

150.465 

325.3213 

<0eXe 

2.504 

1.63608 

0.3842 

0.3829 

1.07742 

Be 

0.059439 

0.059589 

0.042515 

0.041720 

0.082107 

«e 

5.604E-4 

4.891E-4 

1.507E-4 

1.41E-4 

3.187E-4 

De 

3.079E-8 

3.013E-8 

1.351E-8 

1.28E-8 

2.092E-8 

2.6776 


3.20 


2.28107 


Dis. 

En.  2835 


19663.16  35550 


34704 


16055.981 


TABLE  13 


Data  for  Molecular  Bromine 
Part  II.  (79,81)  and  (81,81) 


STATE 

(79,81)B  (2) 

( 79 , 81 ) X  (2) 

(81 ,81)B  (2) 

(81 ,81 )X  (2) 

Te 

15902.47 

0 

15902.47 

0 

“e 

166.5688 

323.3069 

165.5244 

321.29 

“>exe 

1.6159 

1.0641 

1.5957 

1.064 

Be 

.058853 

.081093 

.058118 

.080088 

“e 

4.8007  E-4 

3.1285  E-4 

4.7109  E-4 

3.19  E-4 

De 

2.939  E-8 

2.041  E-8 

2.866  E-8 

1.99  E-8 

Re 

2.6776 

2.28107 

2.6776 

2.28107 

Dis . 
En . 

•* 

19663.7 

* 

16056.988 

* 

19664.184 

** 

16057.996 

All  values  in  inverse  cm  except  Re  which  is  in  Angstroms 
Numbers  in  parentheses  represent  bibliographic  entry  of 
source  article 
*  See  Equation  22 

**  Extrapolated  from  available  data  (see  Appendix  C) 


The  Franck-Condon  Factors  for  molecular  bromine  are 
listed  in  Tables  14  through  19.  Literature  values  of  the 
Franck-Condon  Factors  are  available  for  two  transitions,  E  to 
B  and  D'  to  A'.  The  calculated  and  literature  values  of  the 
Franck-Condon  Factors  for  the  D’  to  A'  transition  are  shown 
in  Table  20.  The  entries  are  in  the  following  order: 
literature  value,  calculated  value,  and  the  absolute 
difference  between  the  two  values.  Although  the 
Franck-Condon  Factors  differ  between  the  calculated  and 
literature  values,  for  the  most  part  the  agreement  between 
both  sets  of  numbers  is  good.  The  most  noticable  derivation 
occurs  at  the  values  of  the  higher  quantum  numbers,  where  the 
IPA  routine  would  have  the  most  effect.  The  potentials  for 
the  literature  values  were  calculated  with  only  the  RKR 
approximation  and  thus  do  not  have  the  added  accuracy  the  IPA 
adds  to  the  higher  energy  region  of  the  curve.  It  is 
believed  that  the  values  presented  in  this  paper  are  more 
accurate  in  this  region  because  of  the  addition  of  the  IPA. 
The  IPA  increases  the  accuracy  of  the  potential  energy  values 
by  repeated  comparisons  between  the  calculated  and 
spectroscopic  energy  values  and  associated  corrections  to  the 
potential  energy  function.  The  most  important  thing  to  note 
from  the  comparison  between  the  calculated  and  literature 
factors  however,  is  that  the  Condon  Parabola  shows  good 
agreement  between  the  two  sets  of  Franck-Condon  Factors. 

Both  sets  of  factors  therefore  predict  the  same  transitions 
to  be  most  probable  even  though  the  factors  themselves  differ 


TABLE  14 


Franck-Condon  Factors  for  Br2  (79)  Lg  B  to  X  Transitions 


V" 

0 

1 

2 

3 

4 

5 

6 

7 

V* 

0 

000 

000 

000 

000 

000 

000 

000 

001 

1 

000 

000 

000 

000 

000 

000 

001 

004 

2 

000 

000 

000 

000 

000 

002 

005 

012 

3 

000 

000 

000 

000 

001 

005 

012 

027 

4 

000 

000 

000 

001 

003 

010 

024 

043 

5 

000 

000 

000 

002 

007 

019 

037 

052 

6 

000 

000 

001 

004 

014 

029 

046 

049 

7 

000 

000 

002 

008 

021 

038 

047 

035 

8 

000 

001 

004 

013 

019 

043 

040 

016 

9 

000 

001 

007 

019 

036 

042 

026 

003 

10 

000 

002 

010 

026 

041 

036 

012 

000 

11 

000 

004 

015 

033 

041 

025 

002 

006 

12 

001 

006 

021 

039 

037 

013 

000 

016 

13 

001 

009 

028 

041 

028 

004 

005 

022 

14 

002 

013 

034 

040 

018 

000 

013 

022 

15 

003 

017 

038 

036 

009 

002 

020 

015 

Values 

in 

the  table  are 

(FCF 

x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


8 


9 


S 


3 


5 


V* 

0 

002 

004 

009 

017 

029 

045 

064 

084 

1 

010 

020 

035 

051 

065 

073 

071 

056 

2 

026 

044 

059 

064 

055 

036 

013 

000 

3 

046 

060 

057 

038 

014 

001 

006 

026 

4 

057 

052 

029 

006 

001 

016 

034 

036 

5 

050 

027 

004 

003 

021 

034 

026 

006 

6 

022 

033 

004 

7 

009 

001 

019 

032 

018 

001 

006 

024 

8 

000 

013 

030 

020 

002 

006 

023 

021 

9 

005 

025 

024 

004 

004 

021 

020 

003 

10 

017 

027 

009 

001 

.  017 

021 

005 

003 

11 

024 

017 

000 

Oil 

022 

008 

001 

015 

12 

023 

005 

004 

020 

013 

000 

Oil 

019 

13 

014 

000 

014 

018 

002 

005 

018 

008 

14 

004 

004 

019 

008 

001 

015 

013 

000 

15 

000 

013 

016 

001 

009 

016 

003 

005 

Values 

in  i 

the  table  are 

(FCF 

x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 
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•*  v  ■ 


-»>: 


,s 


)  095 

166 

219 

215 

151 

071 

002 

000 

000 

000 

120 

087 

017 

013 

121 

217 

074 

Oil 

000 

001 

• 

l  054 

002 

030 

085 

041 

005 

242 

150 

028 

000 

W  - 

025  064  019  012  072  037  229  226  052 


053  020  008  054  013  055  001  176  284 


030  001  041  017  012  000  042  006  119 


002  025  028  002  039  031  021  013  015 


016  007  034  002 


001  025  016  007 


005  028  001  024 


019  015  006  023 


025  002  019  007 


019  002  021  000 


026  015  028  004  039  000 


028  000  001  029  002  035 


007  016  012  016  011  016 


001  024  024  000  025  000 


013  011  Oil  011  009  015 


021  000  000  020  000  018 


007  012  012  007  013  006  008  009  012  004 


000  019  002  017  002  016  017  000  016  002 


003  017  001  016  001  015  012  006  007  012 


011  008  009  008  009  005  002  014  000  014 


017  001  015  001  015  000  001  011  005  005 


016  001  015  001  012  003  008  003  012  000 


010  006  008  007  005  010  012  000  010  005 


004  012  002  013  000  013  010  005  003  010 


000  014  000  012  002  009  003  010  000  009 


001  011  003  008  007  003  000  010  003  003 


006  006  009  002  011  000  003  005  008  000 


011  001  012  000  011  002  008  000  009  002 


TABLE  17  (cont.) 

Franck-Condon  Factors  for  Br2  (79)  E  to  B  Transitions 


4 

080 

001 

014 

001 

001 

000 

000 

000 

000 

000 

5 

321 

105 

004 

027 

000 

004 

000 

001 

000 

000 

6 

076 

336 

121 

Oil 

044 

000 

009 

000 

002 

000 

7 

013 

049 

332 

123 

026 

062 

002 

016 

001 

003 

8 

010 

005 

035 

310 

108 

053 

073 

Oil 

022 

006 

9 

018 

032 

000 

030 

270 

077 

092 

070 

033 

021 

10 

022 

005 

053 

008 

029 

213 

039 

135 

047 

072 

11 

005 

016 

000 

069 

032 

030 

143 

009 

163 

015 

12 

003 

015 

006 

000 

081 

066 

027 

071 

000 

153 

13 

016 

000 

021 

000 

000 

091 

100 

018 

017 

010 

14 

012 

007 

004 

021 

003 

000 

100 

120 

005 

000 

15 

001 

015 

001 

008 

017 

015 

003 

103 

112 

001 

16 

003 

007 

Oil 

001 

009 

014 

036 

012 

094 

076 

17 

011 

000 

012 

005 

005 

006 

012 

065 

024 

067 

18 

Oil 

004 

004 

Oil 

001 

Oil 

001 

014 

096 

034 

19 

003 

Oil 

000 

008 

006 

000 

018 

000 

019 

120 

20 

000 

009 

006 

001 

010 

002 

001 

024 

008 

027 

21 

004 

002 

010 

001 

005 

009 

000 

001 

032 

025 

22 

009 

000 

007 

006 

000 

007 

006 

002 

000 

042 

23 

008 

004 

001 

009 

002 

002 

007 

003 

007 

000 

24 

003 

008 

000 

005 

007 

000 

005 

005 

002 

018 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson’s  Rule  Approximation 


TABLE  18 


ranck-i 

Condon 

Factors  for 

Br  2 

(81)  B 

to  X 

Transitions 

V" 

0 

1 

2 

3 

4 

5 

6 

7 

V' 

0 

000 

000 

000 

000 

000 

000 

000 

001 

1 

000 

000 

000 

000 

000 

000 

001 

003 

2 

000 

000 

000 

000 

000 

001 

004 

Oil 

3 

000 

000 

000 

000 

001 

004 

Oil 

025 

4 

000 

000 

000 

001 

003 

009 

022 

041 

5 

000 

000 

000 

002 

007 

018 

035 

051 

6 

000 

000 

001 

004 

012 

027 

044 

050 

7 

000 

000 

002 

007 

019 

036 

047 

037 

8 

000 

001 

003 

012 

027 

042 

041 

019 

9 

000 

001 

006 

018 

035 

042 

029 

005 

10 

000 

002 

009 

025 

040 

037 

015 

000 

11 

000 

003 

014 

032 

041 

027 

004 

004 

12 

001 

006 

020 

038 

038 

015 

000 

014 

13 

001 

008 

026 

041 

031 

005 

003 

021 

14 

002 

012 

032 

041 

021 

000 

Oil 

022 

15 

003 

016 

037 

037 

011 

001 

018 

017 

Values 

Values 

in  the  table  are 
calculated  using 

( FCF  x  1000) 

Simpson's  Rule  Approximation 

TABLE  18  (cont.) 

Franck-Condon  Factors  for  Br2  (81)  B  to  X  Transitions 


V" 

8 

9 

10 

11 

12 

13 

14 

15 

V* 

0 

002 

004 

009 

016 

028 

044 

063 

084 

1 

009 

018 

033 

049 

064 

074 

072 

057 

2 

024 

042 

058 

065 

057 

038 

015 

001 

3 

044 

059 

059 

042 

017 

001 

005 

024 

4 

056 

054 

033 

008 

000 

014 

033 

037 

5 

051 

030 

006 

002 

019 

034 

028 

008 

6 

032 

007 

002 

020 

033 

022 

003 

004 

7 

011 

000 

017 

032 

020 

002 

005 

023 

8 

000 

010 

029 

022 

003 

005 

022 

022 

9 

003 

023 

026 

006 

002 

020 

0021 

004 

10 

014 

027 

Oil 

000 

015 

022 

006 

002 

11 

023 

019 

001 

009 

022 

010 

000 

014 

12 

024 

007 

002 

019 

015 

000 

009 

019 

13 

016 

000 

012 

019 

003 

004 

018 

009 

14 

006 

003 

019 

010 

000 

014 

014 

000 

15 

000 

011 

017 

001 

007 

017 

004 

004 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


TABLE  19 


Franck-Condon 

Factors 

for 

Br2  (79 

,81) 

B  to  X 

Transitions 

V" 

0 

1 

2 

3 

4 

5 

6 

7 

V' 

0 

000 

000 

000 

000 

000 

000 

000 

001 

1 

000 

000 

000 

000 

000 

000 

001 

004 

2 

000 

000 

000 

000 

000 

001 

004 

012 

3 

000 

000 

000 

000 

001 

004 

012 

026 

4 

000 

coo 

000 

001 

003 

010 

023 

042 

5 

000 

000 

000 

002 

007 

018 

036 

051 

6 

000 

000 

001 

004 

013 

028 

045 

049 

7 

000 

000 

002 

008 

020 

037 

047 

036 

8 

000 

001 

004 

012 

028 

042 

040 

018 

9 

000 

001 

006 

018 

035 

042 

028 

004 

10 

000 

002 

010 

025 

040 

036 

014 

000 

11 

000 

004 

015 

032 

041 

026 

003 

005 

12 

001 

006 

021 

038 

038 

014 

000 

015 

13 

001 

009 

027 

041 

030 

005 

004 

022 

14 

002 

012 

033 

041 

020 

000 

012 

022 

15 

003 

017 

038 

036 

010 

002 

019 

016 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


TABLE  20 


Franck-Condon  Factors  for  Br2 ,  D'  to  A'  Transitions 
Comparison  of  Literature  and  Calculated  Values 

(FCF  x  1000) 


V" 

0 

1 

2 

3 

4 

5 

6 

7 

vr 

- 

- 

0 

0 

0 

2 

12 

44 

0 

0 

0 

0 

0 

0 

2 

10 

38 

- 

- 

0 

0 

0 

0 

2 

6 

- 

- 

0 

0 

2 

11 

39 

94 

1 

0 

0 

0 

0 

2 

10 

37 

88 

- 

- 

0 

0 

0 

1 

2 

6 

— 

— 

0 

1 

6 

26 

69 

H2il 

2 

0 

0 

0 

1 

6 

25 

66 

Cloo) 

_ 

0 

0 

0 

1 

3 

3 

— 

— 

0 

3 

14 

46 

84  | 

67 

3  0 

0 

0 

3 

14 

44 

68 

_ 

_ . 

0 

0 

0 

2 

4 

1 

— 

'  — 

0 

6 

26 

63 

1  73  1 

22 

4 

0 

0 

1 

6 

25 

60 

25 

- 

- 

1 

0 

1 

3 

i 

3 

— 

_ 

2 

12 

39 

|69  | 

46 

0 

5 

0 

0 

2 

11 

38 

47 

1 

_ , 

0 

1 

1 

2 

1 

1 

Order  of  entries:  literature  Franck-Condon  Factor 

calculated  Franck-Condon  Factor 
absolute  difference  of  the  two  values 


75 


TABLE  20  (cont.) 


• 

Franck- 

■Condon  Factors  for 

Br„  , 

D*  to 

A'  Transitions 

Comparison 

of  Literature  and 

Calculated 

Values 

(FCF  x 

1000) 

V" 

8 

9 

10 

11 

12 

13 

14 

15 

• 

V' 

118  ! 

225 

1  f286l  [ 

22lJ 

83 

7 

1 

0 

0 

103 

120 

28 

1 

1 

15 

27 

21 

12 

37 

21 

0 

1 

• 

133 

79 

0 

117 

1  302  ! 

196 

17 

6 

1 

1 

< 

no) 

91 

6 

57 

234 

<5*9) 

86 

2 

3 

12 

6 

60 

68 

53 

69 

4 

c 

61 

0 

58 

64 

9 

266  | 

285 

17 

2 

67 

3 

41 

80 

3 

114  ( 

f309) 

160 

• 

6 

3 

17 

16 

6 

152 

24 

143 

5 

31 

60 

0 

53 

2 

217 

|  330  | 

3 

9 

21 

64 

8 

33 

38 

30 

<T304) 

• 

4 

10 

4 

8 

20 

36 

187 

26 

8 

55 

10 

31 

33 

10 

6 

196 

c 

4 

4 

50 

20 

13 

49 

0 

33 

4 

4 

5 

10 

18 

16 

10 

27 

192 

in  value.  The  coincidence  of  the  parabolas  indicates  that 
the  calculation  of  potential  energy  curves  by  the  RKR  method 
is  adequate  for  most  purposes.  In  Table  20,  the  parabola  is 
represented  by  lines  connecting  maximum  values  within  one 
set.  Thus,  even  though  the  calculated  values  differ  from  the 
ones  published  in  the  literature,  the  transitions  displaying 
maximum  probabilities  have  not  been  shifted  but  show  good 
agreement  with  the  published  maximums.  Similiar  results  are 
seen  when  the  two  sets  of  Franck-Condon  Factors  for  the  E  to 
B  transition  are  compared  and  the  Condon  Parabolas  drawn  in. 
Each  of  the  three  B  to  X  transitions  also  exhibit  clear 
Condon  Parabolas. 

Further  computer  runs  were  made  to  test  the  sensitivity 
of  the  wavefunction  values  to  the  number  of  points  at  which 
they  are  evaluated.  The  potentials  and  waveunctions  for  the 
E  to  B  transition  were  recalculated  using  631  points  to 
evaluate  the  wavefunction.  Program  FCFACT  was  rerun  using 
these  values.  No  substantial  change  in  the  values  of  the 
Franck-Condon  Factors  was  noted.  This  observation  indicates 
that  the  initial  use  of  301  points  insures  accuracy  in  the 
calculation  of  the  wavef unctions . 


Bromine  Chloride 

Two  isotopic  configurations  of  BrCl  were  studied,  ^r^l 
and  "Br^l.  The  Franck-Condon  Factors  for  the  B  to  X 


transitions  were  calculated  for  each  molecule  of  the  pair. 


The  molecular  constants  were  der  ved  from  two  sources,  an 


• 

article  by  J.  A.  Coxon  published  in  the  Journal  of  Molecular 
Spectroscopy  (7)  and  a  thesis  prosectus  written  by  Glen 
Perram  (19).  Because  of  the  lack  of  published  data,  the 
dissociation  energy  must  once  again  be  extrapolated  from 
available  data  (see  Appendix  C).  The  data  for  both  molecules 
of  BrCl  is  listed  in  Table  21.  The  potentials  and  values  of 
the  wavef unctions  were  evaluated  at  301  points  on  each 
vibrational  level.  Here  again,  the  number  of  points  at  which 
the  wavef unctions  were  evaluated  was  increased  to  test  its 
effect  on  the  values.  The  number  of  points  was  increased  to 
901,  almost  triple  the  original  values,  but  the  change  in  the 
Franck-Condon  Factors  was  very  small.  The  Franck-Condon 
Factors  for  both  transitions  studied  are  tabulated  in  Tables 
22  and  23. 

When  the  calculated  results  were  compared  to  the 
literature  values  of  the  Franck-Condon  Factors  for  BrCl, 
there  was  poor  agreenent  between  the  two  sets  of  data.  Table 
24  contains  a  comparison  between  literature  and  calculated 
values  similiar  to  those  described  earlier.  As  can  be  seen 
from  this  table,  not  only  is  the  comparison  between  values 
poor  but  the  maximum  value  lines  do  not  correspond.  These 
lines  should  be  the  Condon  Parabola  but  that  feature  is 
poorly  formed  in  these  two  sets  of  Franck-Condon  Factors. 

The  calculations  for  the  constants  for  the  B  state  in  the 
literature  (7)  include  a  perturbation  term  which  was  not 
included  as  part  of  the  data  for  the  calculations  of  the  B 
state  potential  and  wavef unctions  for  this  thesis.  It  was 


TABLE  21 


Data  for  Bromine  Chloride 

STATE  (79,35)B  (7)  (79,35)X  (7)  (81,35)B  (7)  (81,35)X  (7) 


16879.9 


222.68 


2.88 


107704 


2.5415 


0 


444.276 


3 


.1524695 


16881.2 


221.87 


2.848 


.106919 


2.] 

138 

2.5411 


18945.2 


0 


442.589 


.879 


.1513093 


1.0078  E-7  7.183  E-8  9.9318  E-8  7.074  E-8 


38 


18211.844 


All  values  in  inverse  cm  except  Re  which  is  in  Angstroms 
Numbers  in  parentheses  represent  bibliographic  entry  of 
source  article 

**  Extrapolated  from  available  data  (see  Appendix  C) 


TABLE  23 


anck- 

•Condon 

Factors 

i  for 

BrCl 

(81)  B 

to  X 

Transitions 

V" 

0 

1 

2 

3 

4 

5 

6 

7 

V’ 

0 

000 

000 

000 

000 

000 

000 

001 

001 

1 

000 

000 

000 

000 

010 

2 

000 

000 

001 

002 

005 

010 

020 

033 

3 

000 

001 

003 

007 

016 

030 

046 

059 

4 

000 

003 

009 

021 

038 

054 

062 

054 

5 

002 

009 

024 

044 

060 

061 

044 

018 

6 

005 

022 

046 

064 

060 

036 

009 

000 

7 

014 

043 

067 

063 

033 

005 

002 

021 

8 

030 

068 

071 

037 

004 

004 

025 

032 

Values  in  the  table  are  (FCF  x  1000) 

tr  *  i  “i i  -  i. i  j c  j _ _ f  _  n ..  i  _  i _ _ i 
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not  possible  to  avoid  the  effects  of  this  perturbation  term 
with  the  present  form  of  the  potential  modelling  program. 
Because  of  this,  further  calculations  of  Franck-Condon 
Factors  for  BrCl  were  discontinued. 

Bromine  Monofluoride 

BrF  is  one  of  the  least  studied  molecules  of  the 
interhalogens.  Because  of  this,  the  amount  of  data  available 
in  the  literature  is  small.  Only  one  transition,  the  B  to  X 
has  been  studied.  Values  for  the  molecular  constants  have 
been  reported  for  the  B  state  with  a  natural  mixture  of  the 
^Br  and  fl1Br  isotopes.  The  lower  state  constants  have  been 
reported  for  the  isotopically  pure  ^Br-1^.  These  constants 
were  used  in  the  present  study  to  calculate  the  potentials 
for  the  B  and  X  states  and  Franck-Condon  Factors  for  the  B  to 
X  transition.  The  values  for  these  constants  are  shown  in 
Table  25.  The  wavefunction  was  evaluated  at  301  points 
between  the  minimum  and  maximum  internuclear  distances.  The 
Frank-Condon  Factors  for  the  BrF  B  to  X  transition  are 
included  as  Table  26.  No  literature  values  for  the 
Franck-Condon  Factors  were  uncovered  during  the  literature 
search  conducted  for  this  thesis. 

Iodine  Chloride 

Although  many  states  of  iodine  chloride  are  known,  few 
of  them  have  been  characterized  in  depth.  The  literature 
search  for  this  thesis  produced  molecular  constants  for  five 


18650 


0 


CO 

e 

365 

669.823 

CD  X 

6 

3.753 

e  e 

B 

e 

0.25659 

0.3558193 

a 

2.5826  E-3 

e 

R 

2.0689 

1.758995 

e 

Dissociation 

Energy 


5000 


20622 


4 


TABLE  26 

Franck-Condon  Factors  for  BrF  Large  B  to  X  Transitions 


V” 

0 

1 

2 

3 

4 

5 

6 

7 

V' 

0 

000 

000 

001 

002 

004 

007 

013 

020 

1 

000 

002 

006 

013 

023 

035 

049 

059 

2 

003 

010 

024 

041 

056 

065 

064 

052 

3 

Oil 

032 

056 

070 

067 

049 

025 

006 

4 

032 

068 

081 

064 

032 

007 

000 

Oil 

5 

074 

097 

066 

021 

000 

009 

027 

034 

6 

129 

088 

021 

000 

020 

036 

029 

010 

7 

176 

.  040 

000 

027 

040 

021 

002 

004 

8 

191 

002 

028 

044 

017 

000 

012 

025 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 
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TABLE  26  (cont.) 

Franck-Condon  Factors  for  BrF  Large  B  to  X  Transitions 


V" 

8 

9 

10 

11 

12 

13 

14 

15 

V’ 

0 

030 

041 

052 

064 

074 

081 

085 

085 

1 

066 

065 

058 

045 

029 

014 

004 

000 

2 

034 

016 

003 

000 

007 

020 

033 

040 

3 

000 

007 

021 

033 

035 

027 

013 

003 

4 

027 

034 

027 

013 

002 

001 

010 

022 

5 

024 

008 

000 

005 

018 

025 

022 

010 

6 

000 

007 

020 

024 

016 

004 

000 

007 

7 

018 

024 

015 

003 

001 

010 

019 

018 

8 

019 

005 

000 

010 

019 

016 

005 

000 

Values 

Values 

in  the  table  are 
calculated  using 

(FCF 

Simps 

x  1000) 

on's  Rule  Approximation 
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well-known  states  of  IC1  split  between  two  isotopic 
configurations.  Three  electronic  states,  A,  B,  and  X,  were 
modelled  for  1^21  and  two  states,  B  and  X,  were  modelled  for 
I’fcl.  The  molecular  constants  used  to  calculate  the 
potential  energy  curves  are  shown  in  Table  27.  Not  all  of 
the  constants  for  one  electronic  state  were  taken  from  the 
same  journal  article,  where  available,  the  most  recent 
constants  were  used  in  the  data.  This  effects  primarily  the 
ground  state  of  1^1  because  the  journal  article  from  which 
the  data  was  taken  (the  more  recent  article)  uses  the  older 
constants  as  a  part  of  its  data  set.  The  disssociation 
energy  for  the  B  states  of  both  isotopes  was  calculated  using 
energy  values  averaged  from  both  molecules.  The  authors 
state  their  doubts  about  the  validity  of  this  number  and 
believe  its  true  value  is  within  five  standard  deviations 
rather  than  one  as  normally  quoted  (14:1304). 

The  Franck-Condon  Factors  for  the  three  transitions 
considered  follow  in  Tables  28  through  30.  Each  set  of 
Franck-Condon  Factors  exhibits  the  characteristic  Condon 
Parabola,  although  the  branches  of  the  parabola  are  fairly 
far  apart.  This  represents  a  placement  of  the  states  such 
that  the  equilibrium  internuclear  distances  are  far  apart. 
Literature  values  of  the  Franck-Condon  Factors  are  available 
for  only  a  small  portion  of  one  transition,  the  B  to  X 
transition  of  the  (127,33)  isotope.  The  comparison  of 
literature  and  calculated  values  is  shown  in  Table  31.  As 
before,  the  entries  in  order  of  occurance  are,  the  literature 


TABLE  27 

Data  for  Iodine  Chloride 


STATE 

(127,35) 

A  (8) 

(127,35) 

B  (14) 

(127,35) 

X  (8) 

(127,37) 

B  (14) 

(127,37) 

X  (14) 

T 

e 

13742.9 

17375.58 

0 

17375.47 

0 

0) 

e 

211.0 

204.271 

(14) 

384.323 

200.169 

376.093 

0)  X 
e  e 

2.12 

2.4932 

(14) 

1.5172 

2.5762 

1.4522 

B 

8 

0.08529 

0.086523 

0.114157 

0.082809 

0.109337 

a 

8 

7.4  E-4 

1.402E-3 

5.32  E-4 

1.077E-3 

5.369E-4 

R 

8 

2.685 

2.6659 

2.32091 

2.6667 

— 

Dis . 
En. 

3814.7 

1064.34 

17557.6 

1064.34 

17557.511 

All  values  in  inverse  cm  except  which  is  in  Angstroms 

Numbers  in  parentheses  represent  bibliographic  entry  of 
source  article  for  the  number 


2  OOP  OOP  001  005  016  036  062  080  071  036  005  003 

3  OOP  001  004  013  033  057  071  056  019  OOP  019  047 

4  OOP  002  009  026  050  065  051  015  001  025  047  031 

5  OOP  004  016  040  060  052  019  OOP  023  043  023  OOP 

6  001  007  027  052  056  028  001  015  039  022  OOP  018 

7  002  013  038  056  040  006  006  033  027  001  015  033 

8  003  020  047  052  020  OOP  022  032  005  008  031  015 

9  005  028  053  040  005  008  032  016  001  025  020  OOP 

10  008  036  053  024  OOP  021  027  002  013  026  003  011 

11  012  045  048  010  005  029  013  002  024  011  003  024 

12  018  052  038  001  016  027  002  013  021  OOP  016  018 

13  024  055  025  001  025  017  001  022  009  004  022  004 

14  032  055  013  007  028  005  008  021  001  015  013  001 

15  040  050  004  016  023  OOP  018  012  003  019  002  010  | 

16  049  042  OOP  024  014  003  021  002  011  013  001  017 

17  057  031  002  028  005  010  017  OOP  017  004  008  014 

18  064  020  009  025  OOP  018  008  005  015  OOP  015  005 

19  069  010  017  018  002  020  001  013  008  004  014  OOP 

20  074  003  024  010  008  016  OOP  016  001  011  007  004 

21  076  OOP  029  003  015  009  005  014  OOP  014  001  010 

22  075  002  028  OOP  019  002  011  007  005  011  001  013 

23  072  008  023  002  018  OOP  015  001  011  004  006  009 

24  066  018  015  008  013  002  014  OOP  013  OOP  011  003 


0  126 


1  013 


2  028 


3  048 


4  004 


5  015 


12 


128  118  099  049 


000  012  043  100 


056  060  038  001 


022  001  011  056 


006  033  044  001 


038  028  003  035 


030  017  009 


102  087  064 


021  056  084 


040  011  000 


012  042  052 


038  013 


004  002  001  000 


042  025  013  006 


094  084  064  042 


020  055  082  089 


031  005  003  030 


ES 


16 

017 

EEE 

019  017 

EH5J 


034  009 


GEa 

019 

032 

Oil 

001 

024  « 

023 

023 

002 

010 

031 

018  1 

7  013  001  023  030  005  029  028  004 


022  026  003  030  016 


9  017  025  003  008  011  001  023  023  002  010  031  018 


10  026  006  006  025  002  023  017  000  017  026  005  005 


11  011  002  022  013  021  016  000  018  021  001  011  027 


016  017  000  017  000  016  018  000  013  022  003 


13  008  020  002  010  001  012  018  001  012  020  001  010 


K'EI 


14 

018 

008 

003 

019 

007 

019 

002 

EBEI 

019 

002 

010 

020 

15 

016 

015 

on 

018 

006 

005 

018 

003 

007 

019 

002 

16  005  006  016  000  012  001  01 


17  000  015  006  004  001  011  Oil 


18  006  014  000  013  004  015  001  010  011  000  015  0 


19  013  005  005  013  013  006  004  014  001  010  010  001 


20 

013 

EjI'J 

012 

005 

012 

012 

005 

004 

013 

Esa 

010  j 

21 

006 

003 

012 

000 

005 

on 

EHJ 

012 

004 

005 

012  1 

22 

010 

005 

004 

IsMJ 

012 

003 

006 

010 

K'l'J 

012 

003 

23  001 


24  007 


011  000  010  005  010  000  011  002  007  008  001 


007  002  010  010  003  006  008  001  011  001  008 


TABLE  29 


Franck-Condon  Factors  for  IC1  (35)  B  to  X  Transitions 


V" 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

V' 

0 

000 

000 

000 

000 

002 

006 

015 

033 

061 

091 

114 

127 

1 

000 

000 

000 

002 

007 

019 

043 

071 

091 

085 

055 

020 

2 

000 

000 

001 

005 

015 

037 

062 

075 

057 

020 

000 

013 

3 

000 

000 

002 

009 

025 

049 

063 

047 

013 

001 

025 

050 

4 

000 

000 

003 

013 

033 

053 

048 

017 

000 

023 

045 

031 

5 

000 

001 

004 

017 

038 

048 

028 

002 

012 

039 

029 

003 

6 

000 

001 

006 

020 

038 

038 

012 

002 

028 

034 

007 

004 

7 

000 

001 

007 

021 

036 

027 

003 

00? 

033 

018 

000 

021 

8 

000 

001 

007 

021 

032 

018 

000 

017 

029 

005 

008  029 

9 

000 

001 

007 

020 

026 

Oil 

001 

022 

020 

000 

018 

024 

10 

000 

001 

007 

019 

022 

006 

003 

023 

012 

002 

024 

014 

11 

000 

001 

007 

017 

017 

003 

006 

021 

006 

006 

024 

005 

12 

000 

001 

007 

016 

014 

001 

008 

018 

002 

Oil 

020 

001 

13 

000 

001 

007 

014 

011 

000 

009 

015 

000 

014 

014 

000 

14 

000 

001 

006 

012 

008 

000 

009 

012 

000 

015 

009 

002 

15 

000 

001 

006 

Oil 

007 

000 

009 

009 

001 

015 

005 

005 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


TABLE  29  (cont.) 

Franck-Condon  Factors  for  IC1  (35)  B  to  X  Transitions 


V"  12  13  14  15  16  17  18  19  20  21  22  23 


V' 

0 

128 

118 

100 

077 

034 

020 

Oil 

006 

003 

001 

000 

000 

1 

001 

006 

032 

067 

104 

096 

077 

054 

035 

021 

Oil 

006 

2 

044 

062 

053 

025 

005 

031 

065 

088 

093 

083 

064 

044 

3 

044 

015 

000 

016 

056 

038 

010 

000 

017 

049 

077 

089 

4 

004 

005 

032 

045 

002 

008 

036 

053 

041 

014 

000 

012 

5 

008 

034 

035 

009 

026 

042 

027 

003 

006 

033 

050 

039 

6 

029 

030 

006 

005 

033 

009 

002 

025 

040 

024 

002 

007 

7 

030 

008 

003 

027 

004 

006 

030 

029 

006 

003 

027 

038 

8 

014 

000 

019 

027 

005 

028 

022 

001 

009 

031 

024 

002 

9 

002 

009 

026 

010 

023 

022 

001 

010 

028 

016 

000 

015 

10 

001 

020 

018 

000 

024 

004 

006 

025 

014 

000 

016 

028 

11 

006 

023 

007 

004 

Oil 

001 

020 

018 

000 

013 

024 

007 

12 

014 

018 

000 

012 

001 

010 

021 

003 

007 

022 

009 

001 

13 

018 

Oil 

001 

018 

001 

018 

Oil 

001 

018 

014 

000 

014 

14 

019 

005 

006 

018 

007 

018 

002 

008 

019 

003 

006 

020 

15 

017 

001 

Oil 

014 

013 

Oil 

000 

015 

Oil 

000 

016 

012 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


TABLE  30 


Franck- 

■Condon 

Factors  for 

IC1 

(37)  B 

to  X 

Transitions 

V" 

0 

1 

2 

3 

4 

5 

6 

7 

V’ 

0 

000 

000 

000 

000 

002 

006 

014 

031 

1 

000 

000 

001 

003 

009 

023 

048 

077 

2 

000 

000 

002 

007 

021 

044 

069 

077 

3 

000 

001 

004 

015 

035 

057 

063 

040 

4 

000 

002 

008 

025 

047 

057 

039 

008 

5 

000 

003 

015 

037 

054 

044 

013 

001 

6 

001 

006 

023 

048 

053 

025 

000 

015 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 
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TABLE  30  (cont.) 

Franck-Condon  Factors  for  IC1  (37)  B  to  X  Transitions 


V" 

8 

9 

10 

11 

12 

13 

14 

15 

V' 

0 

056 

085 

110 

126 

129 

122 

104 

081 

1 

096 

088 

057 

021 

001 

006 

033 

068 

2 

054 

017 

000 

018 

048 

062 

047 

018 

3 

007 

004 

033 

053 

037 

008 

002 

025 

4 

003 

032 

047 

022 

000 

014 

040 

039 

5 

025 

040 

018 

000 

021 

040 

022 

001 

6 

036 

019 

000 

019 

036 

015 

000 

021 

Values  in  the  table  are  (FCF  z  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


TABLE  31 


Franck-Condon  Factors  for  IC1,  B  to  X  Transitions 
Comparison  of  Literature  and  Calculated  Values 

(FCF  x  1000) 


8  9 


• 

V'  0 

0 

0 

1  2  7 

18 

37 

63 

naiinJLmJ 

0  0 

0 

0 

0  2  6 

16 

36 

65 

(96Xi20)(i3l) 

►  'w 

c 

0 

0 

0 

1  0  1 

2 

1 

2 

4  3  0 

>» 

o’ 

ST 

0 

0 

1 

3 

10 

25 

48 

73 

86 

77 

49 

17 

1  0 

0 

0 

2 

6 

19 

41 

69  (88) 

81 

52 

■’».  •*.  * 

is 

*  %*  ■ «' 

# 

0 

0 

1 

1 

4 

6 

7 

4 

2 

4 

3 

1  s 

0 

0 

2 

8 

22 

44 

64 

dz] 

45 

14 

0 

15  m 

2  0 

0 

1 

4 

15 

35 

61 

CzD 

56 

20 

0 

13  iy 

• 

0 

0 

1 

4 

7 

9 

3 

7 

11 

6 

0 

•  a 

0 

1 

4 

14 

32 

50  UiJ  34 

7 

2 

22 

42 

3 

0 

0 

2 

8 

24 

48  (64)  49 

14 

1 

24 

49 

0 

1 

2 

6 

8 

2  10  15 

7 

1 

2 

7 

Order 

of 

entries : 

literature  Franck-Condon  Factor 
calculated  Franck-Condon  Factor 
absolute  difference  of  the  two 

(5) 

(14) 

values 

value  of  the  Franck-Condon  Factor,  the  calculated  value,  and 
the  absolute  difference  between  the  two.  The  Condon 
Parabolas  outlined  in  Table  31  show  excellent  agreement 
between  the  two  sets  of  factors.  Even  the  agreement  between 
the  values  of  the  Franck-Condon  Factors  themselves  is  very 
good.  The  comparison  is  only  for  low  vibrational  levels 
which  can  account  for  part  of  the  good  agreement.  The 
closeness  of  the  values  is  even  more  significant  when  it  is 
considered  that  the  two  sets  of  factors  were  calculated  from 
different  sets  of  molecular  constants.  The  literature  values 
are  from  a  1976  article  by  Clyne  and  McDermid  (5)  using 
costants  obtained  from  their  own  spectroscopic  observation 
and  search  of  the  open  literature  while  the  calculated  values 
used  molecular  constants  from  a  1982  article  by  Hansen  and 
others  (14).  The  conclusion  drawn  from  this  is  that  the  B  to 
X  transition  of  the  (127,35)  isotope  of  IC1  is  well  known  for 
the  lower  vibratinal  levels  and  that  further  refinement  of 
molecular  constants  for  these  two  states  will  yield 
improvement  primarily  in  the  higher  vibrational  levels.  The 
combination  of  more  recent  molecular  constants  plus  the  IPA 
routine  used  the  calculation  of  the  potentials  for  this 
thesis  should  lead  to  values  of  Franck-Condon  Factors  which 
are  very  good  at  the  lower  levels  and  very  reliable  at  the 
higher  vibrational  levels. 

Iodine  Monofluoride 

Franck-Condon  Factors  for  transitions  between  two 
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» 


ft 


% 


C 


L 
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electronic  states  of  IF  were  calculated  from  values  for  the 
molecular  constants  found  in  the  literature.  The 
wavef unctions  for  the  B  and  X  states  were  calculated  from 
molecular  constants  found  in  an  article  by  Clyne  and  McDermid 
(5).  The  molecular  constants  found  in  the  article  are  a 
compilation  of  values  from  other  journal  articles.  The 
molecular  constants  are  listed  in  Table  32  and  the 
Franck— Condon  Factors  for  the  B  to  X  transitions  are 
tabulated  in  Table  33.  The  values  of  the  wavef unctions  were 
evaluated  at  301  points  between  the  minimum  and  maximum 
internuclear  distances.  Literature  values  of  Franck-Condon 
Factors  are  included  in  the  article  by  Clyne  and  McDermid. 
These  literature  values  were  calculated  using  wavefunction 
values  based  on  purely  RKR  potentials.  A  side  by  side 
comparison  of  the  literature  and  calculated  values  reveals 
good  agreement  between  the  Franck-Condon  Factors  themselves 
and  very  good  agreement  between  the  two  Condon  Parabolas. 

Discussion  of  Results 

The  Franck-Condon  Factors  for  lead  oxide  were  improved 
by  the  use  of  more  accurate  molecular  constants.  Because  the 
previous  set  of  factors  was  calculated  using  values  of  the 
wavef unctions  based  on  potentials  approximated  by  the  RKR-IPA 
program,  the  old  versus  new  comparison  was  very  close.  The 
most  prominent  changes  were  in  the  higher  vibrational  levels 
where  changes  were  expected  due  to  the  improved  constants. 

For  these  comparisons,  it  can  be  assumed  that  any  changes  in 
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TABLE  32 


Data  for  Iodine  Monofluoride 


STATE 

B  (5) 

X  (5) 

T 

19054 

0 

0 

U) 

406.51 

610.51 

e 

00  X 

1.30 

3.10 

*  *  • 

e  e 

.# 

B 

0.2272 

0.2799 

e 

■  ■  * 

a 

0.001398 

0.001887 

♦  v 

e 

• 

R 

2.1189 

1.9089 

o 

*  v  ‘ 

k;-;; 

Dissociation 

Energy 

11544 

22998 

All  values  in  inverse 

cm  except  Rg 

which  is  in  Angstroms 

Numbers  in  parentheses  represent  bibliographic  entry  of 
source  article  for  the  number 
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V  ' 


0 

005 

031 

088 

157 

202 

198 

153 

094 

046 

018 

006 

002 

1 

023 

093 

151 

118 

031 

002 

063 

140 

158 

118 

064 

027 

2 

056 

138 

097 

007 

031 

097 

067 

004 

029 

111 

147 

115 

3 

096 

125 

015 

030 

086 

025 

010 

078 

066 

004 

029 

109 

4 

128 

069 

006 

080 

025 

015 

072 

022 

012 

077 

053 

000 

5 

142 

017 

050 

051 

005 

065 

016 

021 

065 

009 

026 

078 

6 

138 

000 

076 

005 

049 

030 

013 

058 

003 

039 

050 

000 

7 

119 

016 

059 

009 

056 

000 

054 

008 

032 

039 

002 

057 

8 

094 

045 

024 

042 

020 

030 

030 

014 

044 

001 

050 

012 

9 

070 

070 

002 

058 

000 

051 

000 

047 

003 

039 

016 

020 

10 

049 

082 

004 

044 

017 

029 

020 

027 

016 

032 

008 

042 

11 

032 

080 

024 

017 

043 

003 

043 

000 

042 

000 

040 

003 

12 

020 

069 

046 

001 

047 

006 

032 

015 

022 

019 

020 

017 

13 

012 

055 

061 

004 

030 

028 

007 

037 

000 

038 

000 

037 

14 

007 

040 

065 

019 

009 

042 

001 

032 

013 

018 

022 

Oil 

15 

004 

027 

061 

038 

000 

036 

018 

009 

033 

000 

034 

002 

16 

002 

018 

050 

051 

005 

018 

035 

000 

030 

012 

014 

024 

17 

001 

Oil 

039 

055 

019 

004 

036 

Oil 

Oil 

030 

000 

030 

18 

001 

007 

028 

052 

034 

000 

024 

028 

000 

028 

Oil 

Oil 

19 

000 

004 

019 

045 

044 

008 

009 

034 

007 

012 

028 

000 

20 

000 

002 

013 

036 

048 

020 

001 

027 

022 

000 

026 

Oil 

21 

000 

001 

008 

027 

045 

032 

002 

014 

031 

005 

012 

026 

22 

000 

001 

005 

019 

039 

039 

Oil 

003 

028 

018 

001 

025 

23 

000 

000 

003 

013 

032 

042 

021 

000 

017 

027 

003 

012 

24 

000 

000 

002 

009 

024 

039 

030 

005 

006 

027 

014 

001 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 
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TABLE  33  (cont.) 

Franck-Condon  Factors  for  IF  Large  B  to  X  Transitions 


V" 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

V* 

0 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

000 

1 

009 

002 

001 

000 

000 

000 

000 

000 

000 

000 

000 

000 

2 

064 

026 

008 

002 

000 

000 

000 

000 

000 

000 

000 

000 

3 

139 

105 

054 

021 

001 

000 

000 

000 

000 

000 

000 

000 

4 

044 

120 

132 

087 

013 

003 

001 

000 

000 

000 

000 

000 

5 

030 

003 

072 

133 

065 

024 

006 

001 

000 

000 

000 

000 

6 

050 

067 

008 

025 

135 

092 

041 

012 

003 

000 

000 

000 

7 

024 

010 

070 

039 

067 

134 

117 

061 

021 

005 

001 

000 

8 

023 

053 

002 

041 

008 

028 

115 

134 

083 

032 

008 

001 

9 

042 

000 

050 

025 

069 

033 

004 

084 

140 

106 

046 

013 

10 

000 

045 

012 

023 

000 

049 

058 

001 

050 

134 

125 

063 

11 

036 

014 

021 

036 

052 

014 

021 

070 

015 

022 

118 

140 

12 

025 

010 

036 

001 

004 

034 

037 

003 

064 

036 

005 

094 

13 

000 

038 

000 

038 

031 

023 

012 

051 

001 

045 

056 

000 

14 

026 

010 

025 

014 

024 

010 

040 

000 

049 

014 

024 

067 

15 

030 

005 

029 

005 

003 

037 

000 

042 

005 

034 

032 

007 

16 

005 

030 

001 

032 

033 

002 

034 

006 

030 

020 

016 

045 

17 

005 

022 

012 

016 

013 

020 

015 

019 

021 

013 

034 

003 

18 

025 

001 

029 

000 

003 

027 

005 

028 

004 

033 

002 

041 

19 

026 

008 

014 

019 

026 

001 

030 

000 

032 

000 

035 

001 

20 

008 

025 

000 

026 

015 

013 

012 

019 

008 

024 

007 

026 

21 

000 

022 

Oil 

008 

000 

027 

001 

024 

006 

020 

Oil 

019 

22 

Oil 

006 

025 

001 

010 

009 

021 

002 

026 

000 

028 

001 

23 

024 

000 

019 

014 

024 

001 

021 

009 

Oil 

019 

004 

025 

24 

023 

Oil 

004 

024 

015 

016 

002 

024 

001 

021 

007 

014 

Values  in  the  table  are  (FCF  x  1000) 

Values  calculated  using  Simpson's  Rule  Approximation 


the  Franck-Condon  Factors  occurred  because  of  the  new 
constants  and  not  because  of  a  difference  in  the  modelling  of 
the  potential.  The  same  however,  can  not  be  said  of  the 
calculations  done  for  the  interhalogens.  Generally,  the 
potentials  in  the  literature  were  modelled  using  a  straight 
RKR  approximation,  without  the  added  accuracy  of  the  IPA 
routine.  When  the  literature  and  calculated  values  of  the 
turning  points  for  the  potential  are  compared,  the  changes 
induced  by  the  IPA  become  apparent.  The  width  (R(max)  - 
R(min))  of  the  two  energy  levels  are  approximately  the  same 
for  any  vibrational  level  (the  agreement  is  better  at  the 
lower  V  levels)  but  the  values  of  the  turning  points  are 
shifted  towards  smaller  internuclear  distance,  as  can  be  seen 
in  Figure  III.  This  shift  is  not  constant  for  all  energy 
levels  however,  as  might  be  expected  if  one  of  the  rotational 
constants  were  wrong,  but  instead  gets  larger  with  larger 
values  of  V.  This  shifting  of  the  RKR-IPA  potential  from 
straight  RKR  values  for  potentials  calculated  with  the  same 
molecular  constants  shows  the  IPA  routine  at  work.  Because 
the  iterative  processes  of  the  IPA  serve  to  improve  the 
accuracy  of  the  potential  (especially  at  higher  vibrational 
levels)  by  repeated  comparisons  and  corrections  of  the 
potential  energy  values,  these  shifted  potentials  are  more 
accurate  than  those  reported  in  the  literature.  Thus,  the 
wavef unctions  calculated  from  improved  potentials  are  more 
accurate  and  from  them  the  Franck-Condon  Factors  gain 
accuracy.  Also,  more  accurate  potentials  allow  the  number  of 
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FIGURE  III, 


INTERNUCLEAR  SEP.  (Ang) 


Comparison  of  Turning  Points  for  RKR  and 
RKR-IPA  Methods.  Turning  Points  shown  are 
for  BrF  ground  state. 
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V.  Summary  and  Recommendations 

This  thesis  covers  the  modelling  of  potential  energy 
wells  using  a  Rydberg-Klein-Rees  approximation  plus  an 
optional  Inverted  Perturbation  Approach  routine  to  increase 
the  accuracy  of  the  potential  energy  curve  at  higher 
energies.  The  wavefunction  is  approximated  for  each 
rotationless  vibrational  level  by  calculating  its  value  at  a 
specified  number  of  points  along  the  region  of  interest. 
Finally  using  these  tabulated  wavefunction  values,  the 
Franck-Condon  Factors  (which  are  proportional  to  the 
transition  probability)  are  calculated  for  24  transitions 
split  between  a  total  of  six  molecules  (lead  oxide  plus  five 
interhalogens).  The  interhalogens  are  of  special  interest  to 
the  people  at  the  Air  Force  Weapons  Lab  at  Kirtland  AFB. 

This  thesis  opens  with  a  discussion  of  the  RKR 
approximation,  including  when  the  approximation  is  used  and 
the  Dunham  Equation.  This  equation  is  not  only  the  starting 
point  for  calculation  of  the  potential,  but  also  the  equation 
which  defines  the  molecular  constants.  The  discussion 
continues  emphasizing  the  importance  of  the  "f"  and  "g" 
variables  in  the  calculation  of  the  turning  points.  The 
potential  is  modelled  by  a  set  of  these  turning  points 
(maximum  and  minimum  internuclear  spearation  of  the  two 
nuclei)  listed  in  order  of  increasing  energy.  This  ends  the 
first  cut  calculation  of  the  potential.  The  articles 
reviewed  which  published  the  molecular  constants,  stopped 


their  work  here  and  calculated  the  Franck-Condon  Factors 
based  on  an  RER  potential.  The  recently  developed  IPA 
routine  was  used  in  this  thesis  to  increase  accuracy  beyond 
the  RER  analysis.  The  deviation  from  a  straight  RER  approach 
was  most  often  seen  in  the  higher  vibrational  levels  as  a 
slight  shift  of  the  potential  curve  turning  points.  The  IPA 
is  referred  to  as  an  inverted  approach  because  the  AE^  in 
equation  14  is  known  and  AUQ(r)  is  desired,  generally  the 
opposite  is  true.  The  theory  and  the  procedures  of  the  IPA 
follow  the  discussion  of  the  RER  method. 

Next,  the  Franck-Condon  Principle  and  the  calculation  of 
Franck-Condon  Factors  is  studied.  The  principle  is  described 
both  classically  and  quantum  mechanically  and  differences 
between  the  two  are  mentioned.  Also  the  concept  of 
transition  probability  based  on  the  maxima  and  minima  of  the 
wavef unctions  of  the  two  states  is  discussed.  The  formula 
for  the  Franck-Condon  Factors  is  derived  from  the  quantum 
mechanical  statement  of  the  transition  moment  for  two 
electronic  states  of  a  diatomic  molecule.  Once  the  final 
form  of  the  Franck-Condon  Factor  formula  is  obtained  the 
factor’s  relationship  to  the  probability  of  a  transition 
occurring  is  more  easily  seen.  The  section  which  calculates 
the  Franck-Condon  Factors  was  upgraded  from  a  simple 
histogram  approach  to  a  Simpson's  Rule  approximation.  The 
results  of  this  modification  were  verified  by  comparing  the 
newly  calculated  Franck-Condon  Factors  with  a  past  set 
calculated  using  the  same  input  data  and  the  older  version  of 


the  program. 

Finally,  the  Franck-Condon  Factors  themselves  were 
calculated.  Nine  states  were  modelled  for  lead  oxide  and 
factors  for  10  transitions  were  calculated.  For  molecular 
bromine  a  total  of  seven  states  were  modelled  using  three 
isotopic  forms  of  the  molecule.  The  Franck-Condon  Factors 
for  six  transitions  (between  states  of  the  same  isotopic 
make-up)  were  tabulated.  For  bromine  chloride,  two  states 
each  for  two  isotopes  were  approximated  and  the  factors  for 
one  transition  per  molecule  were  calculated.  Bromine 
monofluoride  had  just  two  states  modelled  and  one  set  of 
Franck-Condon  Factors  calculated.  A  total  of  five  electronic 
states  were  approximated  for  iodine  chloride  split  between 
two  isotopic  forms.  Within  these  five  sets  of  potential 
wells  three  transitions  had  the  factors  calculated  for  them. 
Finally,  two  states  and  one  transition  were  calculated  for 
iodine  monofluoride.  With  the  exception  of  bromine  chloride, 
good  agreement  between  the  most  probable  transitions  (the 
Condon  Parabola)  for  each  set  of  factors  was  found  even 
though  the  numerical  values  of  the  two  sets  differed.  This 
shows  that  the  less  accurate  RKR  approach  is  sufficient  for 
many  uses  of  Franck-Condon  Factors. 

Recommendations 

1.  The  rescaling  and  shifting  subroutines  of  program 
FCFACT  need  to  be  modified  and  thoroughly  tested.  It  was 
discovered  late  in  the  quarter  that  these  routines  under 


special  circumstances  do  not  function  as  expected.  The 
problems  are  more  than  likely  minor,  but  the  present  form  of 
the  program  forces  the  two  input  wavef unctions  to  be  scaled 
identically  and  defined  on  the  same  range.  If  this  is  not 
done,  the  program  will  calculate  a  set  of  Franck-Condon 
Factors  which  are  strongly  shifted  towards  the  lower 
vibrational  quantum  numbers.  Allowing  the  two  sets  of 
wavefunction  values  to  be  entered  with  different  scaling 
factors  will  increase  the  power  of  the  program  by  allowing 
the  most  accurate  set  of  wavefunction  values  for  a  particular 
state  to  be  used. 

2.  The  RKR-IPA  program  should  be  modified  to  accept 
data  for  a  potential  which  has  been  calculated  by  other 
means.  The  potential  data  could  be  entered  as  turning  points 
and  energy  values  for  each  vibrational  level  and  could  be 
drawn  from  journal  articles.  The  RKR  portion  of  the  program 
would  be  skipped  and  the  IPA  applied  to  the  already  existing 
potential.  Several  problems  of  working  with  literature  data 
could  be  avoided  this  way.  Examples  of  these  would  be  cases 
where  the  authors  used  more  molecular  constants  than  were 
reported  in  the  article  or  if  the  constants  reported 
contained  an  extra  perturbation  term  which  can  not  be 
duplicated.  This  problem  occurred  for  one  molecule  (BrCl)  in 
this  thesis. 

3.  Matrix  solutions  to  the  Franck-Condon  Factor 
integral  should  be  investigated.  Solutions  using  matrix 
techniques  have  been  proposed  in  the  literature  (12,  1,  25). 


These  methods  would  increase  the  efficiency  of  the 
calculation  and  be  a  good  match  to  programs  currently  being 
developed  which  model  the  potential  wells  using  matrix 
techniques.  These  techniques  become  especially  important  a 
computer  usage  and  demands  for  computer  time  increase. 
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Appendix  A 


FCFACT,  Revised  Code 

This  appendix  contains  the  upgraded  calculation  of  the 
Franck-Condon  Factors  with  the  Simpson's  Rule  approximation 
to  the  overlap  integral.  Changes  to  the  program  included  the 
addition  of  the  file  WFXN  which  holds  the  product  of  the  two 
wavefunction  files.  Other  changes  increased  the  limits  over 
which  the  integral  was  evaluated  to  include  the  entire 
interval  on  which  the  potential  was  modelled.  The  final 
change  made  to  the  program  was  the  upgrade  from  a  simple 
Histogram  approach  to  a  Simpson's  Rule  Approximation  to  the 
overlap  integral.  None  of  the  Scale,  Shift,  or  Normalization 
routines  were  changed.  Only  the  portion  of  the  program  which 
calculates  the  Franck-Condon  Factors  is  included  in  this 
appendix,  the  other  portions  of  the  program  remain  as  listed 
in  Capt  Pow's  thesis  (20). 


a  a 


»ii  :yji] 


BEGIN  CALCULATION  OF  FCF'S.  CONSIDER  ALL  V  LEVEL  CONBINATIQNS.  000877 
FIRST  MULTIPLY  WAVEFUNCTIONS  TOGETHER  FOR  TRANSITION  OF  INTEREST,  000884 


THEN  INTEGRATE  BY  SIMPSON'S  RULE  OVER  ENTIRE  ARRAY.  FINALLY  000891 

CALCULATE  THE  FCF  (SQUARE  OF  THE  INTEGRAL).  000898 

000905 

IF  (NPTLO  .LT.  NPTUP)  THEN  000912 

NPT-NPTLO  000919 

ELSE  000926 

NPT-NPTUP  000933 

ENOIF  000940 

000947 

00  20,  I-O.NVUP  000954 

00  25,  J-O.NVLOW  000961 

DO  30,  K=1 ,NPT  000968 

30  UFXN(K)-)DFXNL0(J,K)*1ilFXNUP(I,K)  000975 

000982 

♦+++  BEGIN  INTEGRATION  -  'am'  HOLDS  VALUE  OF  INTEGRAL  +++♦  000989 

000996 

SUM-(WFXN(1 )  -  liFXN(NPT))/2.  001003 

DO  40,  K-2,(NPT-1),2  001010 

40  SUM-SUM  +  2.*UFXN(K)  +  UFXN(K+1)  001020 

001030 

SUM-2. *SUM/3.  001040 

001060 

♦♦++  CALCULATE  FCF  (SQUARE  OF  THE  INTEGRAL)  001070 

001080 

FACTOR(I,J)-SUM*SUM  001090 

CONTINUE  001100 


Processing  a  Batch  Job 


Programs  are  run  in  batch  mode  either  because  of  the 
large  amount  of  output  generated  or  because  the  program 
requires  large  amounts  of  computer  memory  be  available  to 
run.  The  time  required  to  run  a  batch  job  on  the  Cyber  has 
been  greatly  reduced  with  the  addition  of  one  system  (System 
B)  entirely  dedicated  to  processing  batch  jobs.  If  the  jobs 
are  submitted  during  non-peak  hours  such  as  weekends  or  early 
morning  the  entire  job  can  be  processed  by  the  computer  in 
under  two  minutes.  A  run  during  peak  hours  may  take  as  long 
as  six  hours  to  be  processed. 

Batch  jobs  are  submitted  using  a  submit  command  and  a 
jobfile,  which  is  a  file  containing  the  identification  of  the 
user  and  the  list  of  commands  to  be  processed  by  the  compu 
ter.  The  exact  form  of  the  jobfile  will  vary  from  run  to  run 
and  from  program  to  program.  Two  examples  are  shown  in 
Figures  B-l  and  B-2.  Figure  B-l  shows  the  jobfile  for  the 
RKR-IPA  program  which  must  have  the  input  data  file  read 
external  to  the  main  program.  Figure  B-2  shows  the  jobfile 
for  program  FCFACT.  This  program  contains  internal  calls  for 
the  data  files  and  thus  they  do  not  need  to  be  read  in  by  the 
jobfile.  The  first  four  lines  of  the  jobfile  are  the  same 
for  all  jobfiles  but  the  rest  of  the  file  can  have  several 
variations,  many  of  which  will  accomplish  identical  tasks. 


The  specific  commands  required  and  their  function  are  covered 
in  the  user's  manuals  written  for  the  Cyber  and  its  operating 
system.  Once  the  jobfile  has  been  created  and  saved,  the  job 
is  submitted  using  the  SUBMIT  command,  the  filename  of  the 
jobfile,  and  a  destination  for  the  output  (to  route  output  to 
AFIT,  use  RB-AF).  All  of  the  options  for  this  command  are 
covered  in  the  user's  manual  for  the  Cyber. 


/JOB 

LDB.STCSB. 

USER.T841277 , PASSWD. 
CHARGE,*. 

FTN5.L-0. 

LGO. 

SAVE, WAVE. 

/NOSEQ 

/EOR 

/READ, IPA 
/EOR 

/READ.BRFLGB 

/EOF 


FIG  B-l.  Jobfile  for  RKR-IPA 


/JOB 

LDB.STCSB. 

USER.T841277, PASSWD. 
CHARGE,*. 

GET.BRFBW. 

GET.BRF79XW. 

FTN5.L-0. 

LGO. 

/NOSEQ 

/EOR 

/READ.FCFSMPR 

/EOR 


FIG  B-2 .  Jobfile  for  FCFACT 
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Appendix  C 


Extrapolation  of  Dissociation  Energy  from  Known  Data 


In  several  cases  the  dissociation  energy  was  not 
available  in  the  literature  in  any  form,  either  chemical  or 
spectroscopic.  In  cases  where  the  other  molecular  constants 
were  available  for  more  than  one  molecular  form  of  the  same 
electronic  state  of  the  molecule,  the  values  could  be 
extrapolated.  The  approximation  assumes  similiarity  between 
these  identical  states  of  different  molecular  forms.  To 
calculate  the  desired  dissociation  energy,  the  value  of 
(Dissociation  Energy  +  /2)  was  assumed  to  be  equal  for  both 

forms  (See  equation  below).  The  equation  can  be  set  up  and 
solved  for  the  desired  energy  quickly  and  calculates 
dissociation  energies  consistant  with  known  values. 
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